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1 Introduction 



Thanks to fine work by Chetyrkin [Q, we have information about the imaginary part of 
the correlator of a pair of scalar currents at the 4-Ioop level of QCD, i.e. to the same order 
in perturbation theory as for the vector channel [@]. The role of the vector correlator, in 
the annihilation of an electron-positron pair into hadrons, has been much studied. Here 
we focus on the scalar channel, which is even more intriguing. There are at least 6 potent 
reasons for studying scalar correlators. 

1. Quark masses: Consider the scalar divergence, i{ma—'rnh)'4' a'4'bi of a flavour- changing 
vector current. In perturbative QCD, its correlator is known exactly to two-loop 
order, where it involves trilogarithms of two variables: vn^jQ^ and vn^jQ^ ^ at euclidean 
momentum . At 3 loops, it is was possible to obtain the first two 0, ^ H, terms of 
the expansion in at 4 loops, only the imaginary part is known, and then only to 
leading order in The coefficients of quark and gluon condensates are known, to lesser 
accuracy §]. Combining this information with experimental data on mesons, coupling 
to the scalar form factor of semi-leptonic Kij, decay, and with further theoretical input 
from chiral perturbation theory, the scale-dependent strange-quark mass of the modified 
minimal subtraction (MS) scheme has been determined at a scale of /i = 1 GeV, from 
whence it may be evolved to higher scales, using the anomalous quark-mass dimension 
Imis^s) = d\ogm{fi^)/d\og{^'^). 

2. Higgs decay: At the opposite extreme of very high energy, 4-loop perturbative anal- 
ysis of the scalar channel yields radiative corrections [0 to the decay of the Higgs boson 
of the standard model into quark-antiquark pairs, with a coupling to each fiavour pro- 
portional to the mass of the quark (the 3-loop corrections are known from the results of 
Refs.|, I, 0,0)). 



3. Renormalization group: As is clear from these two important phenomenological 
investigations, mass renormalization is of the essence in the scalar channel. In the vector 
channel, we may ignore quark-mass effects at high energy; in the scalar channel they 
abide, since the only form of scalar current that has meaning is ^/'^Ma^bV'bj where M is 
a mass matrix. Hence the vertex renormalization of the scalar current il^ai^b is precisely 
the inverse of mass renormalization. It follows that the anomalous quark-mass dimension 
7m(cts) is ever present in the renormalization-group equations for correlators of scalar 
currents, while in the vector case it is inactive at very high energy, provided the order 
0{m^ IQ"^) corrections are neglected. 

Thanks to the dedicated labour and great ingenuity of colleagues, we have been pro- 
vided with the 4-loop anomalous quark-mass dimension 0] and the 4-loop beta func- 
tion |T^ of QCD. In our scalar analysis, these are inextricably intertwined. 

4. Renormalons: It is thus of great interest to try to extend our understanding of 
perturbative quantum field theory by studying the interplay of coupling-constant and mass 
renormalization. In this respect, we noticed an apparent omission, concerning behaviour 
at higher orders in the scalar channel. Let Nf be the number of quark fiavours and as(/^^) 
be the strong coupling at scale yU^. In the limit Nf oo, with h = Njas/Qir held fixed. 



1 



the vector correlator is known |T^ to all orders in b, at order 1/Nf. At first sight, this 
limit appears remote from asymptotically free QCD, where the beta function is dominated 
by gluons. However, it has become common practice to transform large- A'^y results to so- 
called large-/3o results, where Pq = (llA'c — 2A^j)/12, with = 3 colours and Nj = 3, 4, 5 
active flavours, gives the one-loop term in the QCD beta function 

dlogas ^ fas 



By the simplistic device Nf — > Nf — called naive nonabelianization (NNA) in |15|, 



one transforms the irrelevant large- A^/ ultraviolet (UV) factorial perturbative growth of 
an abelian theory, like QED, into the highly pertinent large-/3o asymptotic perturbative 
growth of QCD series, since —b is then replaced by 6 = /3oas(/i^)/7r 1/ log(/i^/A^). Then 
the so-called renormalon structure (i.e. the perturbative factorial growth) of the vector 
result [14] is related to the way that long-distance physics is absorbed into condensates in 
the operator-product expansion (OPE) of the vector correlator of QCD (for the studies in 
QCD see e.g. |T|, ^ |18[ and |g for the reviews). 

The virtue of an all-orders large- A^/ result is to provide a smooth map from the full 
two- loop result, which it exactly reproduces, to the large-order behaviour, which it re- 
produces at leading order in 0{1/Nf). It was therefore natural to inquire whether the 
vector analysis [jl4| had yet been extended to the scalar case. We did not find such a 



work. Undertaking the task ourselves, we came to understand why it is so much more 
difficult in the scalar channel: the ultraviolet (UV) infinities of mass renormalization must 
be included to all orders in the coupling. In Sec. 2, we achieve this, after taking guidance 
from the study of critical phenomena. 

5. Multiloop estimates: Effort has been expended in estimating effects in electron- 
positron annihilation, r-decay and in deep-inelastic scattering characteristics beyond the 
orders of perturbation theory that are exactly computed (for different approaches see 
Refs. pT|, E2, pOl). A similar attempt was made in the case of the decay width of the 



Higgs boson |23| following the ideas of Ref.|21] and by the authors of Ref. |24] using the 



asymptotic Pade-approximant method of Ref. P3|. This can only be inspired guesswork, 
informed by past experience and hopeful intuition. To study in detail the ideas and 
approximations, lying beyond such guesswork, one must test them in detail in as many 
processes as possible. In Sec. 3 we submit a variety of suggestions to detailed scrutiny in 
the scalar channel, by taking account of existing 4-loop input from |ll], 12, 13| and by 



exploiting our new all-orders results at large- A^/. Special attention is paid to the results of 
application of the NNA procedure, closely connected to the large- A^j expansion. We also 
investigate the structure of the 3- and 4-loop coefficients we are interested in within the 
"dual NNA" procedure, which is exact in the large number of colours A,, limit. 

6. Analytical continuation: Direct multiloop calculations are usually performed in 
the euclidean space-like region. However, for the scalar correlator the quantities of phe- 
nomenological interest, namely the spectral functions of the QCD sum rules (see e.g. 
Refs. P^, P7| , P^ ) and the decay width of the Higgs boson are proportional to its imagi- 
nary part, defined in the minkowskian time-like region. In higher orders the coefficients of 
perturbative expansions contain vr^-contributions, which generally speaking, are not small 



2 



and can affect the asymptotic structure of the perturbative series in powers of as/vr. These 
effects have been much studied (see e.g. Refs.p9|-p0|,[^]) in attempts to resum them in 



all-orders of perturbation theory. In Sec. 3, using the results of Sec. 2, we develop further 
this approach to the case of fractional powers of a^, which appear in the relation between 
the MS-scheme running quark mass m(/i^) and the scheme-invariant mass m, introduced 
in Refs. PT|, (recent analogous independent considerations were given in Ref . ||43|| ) . The 



NNA approximation and the Borel resummation technique will be essential theoretical 
cornerstones of our analysis. 

In the Conclusions we discuss the theoretical uncertainties of the various approaches 
and results, considered in the previous sections, and summarize the phenomenological 
relevance of the results obtained. 



2 Renormalon analysis at large Nf 

First, we briefly review the much easier vector case, where the Ward identity Zi = Z2 
protects the renormalon chain from UV disturbance, by cancellations between the two- 
loop skeletons into which the chain is inserted. All of these vector methods are necessary 
here, though they are not sufficient. Then we turn to the scalar case, which requires more 
powerful techniques, since there is no Ward identity to protect the insertions. Our aim is 
to handle both UV-divergent two-loop skeletons exactly, and to renormalize the mass to 
all orders in Nfag at large Nf, so as to achieve an MS-renormalized result that connects 
the full two-loop result to the large- A^/ renormalons, checking en route the 0{1/Nf ) terms 
in the 3- loop correlator and 4-loop imaginary part that were given in 

In the interests of transparency, we describe the situation at large- A^/, in the first 
instance, since this is a well-defined limit. Only later do we allow ourselves the prevailing 
luxury of believing that this has anything to do with QCD. By presenting things this way, 
we allow the reader to distinguish hard (and also new) analysis from more easy (and also 
old) conjecture. 



2.1 Combinatorics of vector resummation 

Consider the formal series 

b \ " ^ Lie, ne) 



where the multiloop generating function 

Lie, 6)= L,,ke'S'' (3) 

j,k>0 

is regular near e = S = 0. Here b stands for a renormalized coupling chosen such that the 
large- A^/ beta function vanishes aX b = e in d = 4 — 2e dimensions. In (^), the denominator 
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b — e comes from transforming the bare to the renormahzed couphng. For example, in 
minimally subtracted large- A^/ QED, the bare charge cq is related to b = A^/a(/i^)/37r by 



36 (/iV47r)T(l - 2e) 



47r2 i-b/e [T{l-e)]^T{l+e) 

where /x^ makes b dimensionless, and the F functions make one-loop massless two-point 
diagrams rational in e at euclidean momentum fi^ = Q^. 

Now we resum the series by collecting powers of the renormahzed coupling, obtaining 

n>l k<n ^ 

where the Laurent series at n > 1 loops starts, by assumption, with At any order 

in renormahzed perturbation theory, there are only two significant terms, namely [1^ 

Cn,l = 7T (6) 

n{n — I) 

Cnfl = -^^ + (-ir(n-2)!Lo,„_i (7) 
n[n — 1) 

which are obtained by formal combinatorics. The first result determines the n-loop contri- 
bution to the anomalous dimension; the second gives the finite part. Thus, in this simple 
vector case, protected by the Ward identity Zi = Z2, it is not necessary to know everything 
about the master n-loop integral L{e,ne), resulting from a chain of n — 2 fermion loops 
in a pair of two-loop skeletons. It suffices to know L{e,0) and L{0,6), i.e. the analytic 
continuation to zero loops, and the Borel limit e ^ 0, with 6 = ne fixed. 



2.2 Analysis of vector resummation 



In the case of the correlator of the vector current, the analytical problem was solved in 
closed form by the first author, who obtained [|1^ 



L{e,0) 
L{0,5) 



(l + e)(l-2£)(l-2£/3) 
B{2 - e,2 - e)T{3 - 6)T{1 + e) 



32 



E 



2-6^, (P - (1 - 6W 



(9) 



The simple F-function result @ was long since known from [Q. The all-orders re- 
sult (I) of 0] is in agreement with the independently calculated perturbative expansion 
of Ref.0. At two loops, 

^(0,0) = 16E(^^3Ip =3 (10) 



k=2 



gives the Jost-Luttinger |4^ singularity, when b = a/37r. 



4 



Before explaining how to solve the demanding case of the scalar correlator, we need to 
explain how the far easier vector case was handled. At n-loops, there are two diagrams: in 
the first a chain of n — 2 fermion loops dresses a fermion line; in the second, it is exchanged 
between fermion and antifermion. The first case is easy: only F functions occur; the second 



is hard: there is an F3 2 hypergeometric series of the form given in [0. However, this 
series is needed only in the Borel limit, where it gives a trigamma function. The series 
expansion in 6 of the 4-dimensional two-point two-loop scalar diagram, with a modification 



{1/k of the momentum dependence of the totally internal propagator, is E7 



z>o 



with /(O) = 6(^(3) giving the familiar unmodified result (see e.g. |^). The renormalon 
series in (^) is obtained from the powerful trigamma identity 

which cancels the singularities of (|ll|) at both 6=1 and 6 = 2. The vector result (P) has 
double poles at the negative integers, and at positive integers greater than 2. Yet it has 
no pole at 5 = 1 and merely a single pole at 5 = 2. In the QCD case, the pole at 5 = 2 
signals the need to absorb long-distance effects into a nonperturbative gluon condensate, 
at order l/Q^. No pole can appear at 5 = 1 since there is no dimension-2 gauge-invariant 
operator in the massless theory that could produce nonperturbative effects of order l/Q^- 
The appearance of e^/^ in (^ is easy to understand: in d = 4 — 25 dimensions, each of the 
n — 2 one-loop insertions brings with it a factor 

giving [/(e^)]""^ —>■ exp(55/3), in the Borel limit. To suppress it, one may set yU^ = 
Q^exp(— 5/3), in the MS scheme, corresponding to the perhaps more physical procedure 
of subtracting the one-loop photon propagator at Q^, in the MOM-scheme in QED or 
l^-scheme in QCD (for its 2-loop definition see Ref. 



2.3 Combinatoric and analytic complexity in the scalar case 

In the case of the correlator of the scalar current mipip both the combinatorics and the 
analysis are more demanding. The behaviour at large and large Nf is given by 

n, = U,{e)(^l + 2C,^-^^^±^^ (14) 
U,{e) ^ -2[m{f^')rQ'd,^J^l^^ (15) 
lA{b,e) ^ /f^da: (16) 
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with b = TpN fas^ii^) I'i'n giving the large- iV/ contribution to the beta function. Here Nf 
is the number of fermion flavours, is the value of quadratic Casimir operator in the 
fundamental fermion representation, rfj? is the dimensionality of this representation, and Tp 
specifies the normalization of the corresponding matrices. Concretely, Cp = 4/3, dp = 3, 
and Tp = 1/2, in QCD. In the abelian case of QED, we simply set Cp = dp = Tp = 1. 

The stumbling block is the MS-renormalized mass m^jj?) in the one- loop term (plsl). 
Renormalization of the scalar vertices is mandatory: without it the discontinuity of lis in 
the physical (i.e. minkowski) region —Q^ = s > would be infinite at e = 0. At large Nf, 
the multiplicative renormalization Zm = mo/m(/i^) may be considered as additive, with a 
vertex counterterm Zm — 1 giving integral ([161) . The numerator g of the integrand is the 
all-orders contribution to the anomalous mass dimension 

ctlog(/i^) TpNf ' 

at large N f. The one-loop value (7(0) = 9/4 gives 7m(«s) = —Ois/T^ + 0{a1) in QCD. 

It is straightforward to insert a chain of fermion loops in the one-loop diagram for the 
scalar vertex and obtain the critical exponent 

id -If Tid-2) sinvre , , 



at large Nf m d = A — 2e dimensions, giving the all-orders result for g{b) in (|T^), by virtue 
of the fixed point at b = e, to leading order in 1/Nf. We note, for future reference, that 
g{e) is finite for all > — 1, and hence that g{b) is finite for b < 5/2. The validity of ([T9|), 
for all d > —1, is the true origin of the one-loop result g{0) = 9/A, at d = 4. 

Anomalous dimensions at 0{1/Nf) are given by functions that differ from ( |T9D only by 



multiplication of a rational function of d. The 0(1/ Nj) term in (|TB|) was found in in the 



case of QED, and very recently in ||51| for the yet more demanding nonabelian case of QCD. 
There one finds derivatives of F functions, which still give zeta values. At 0{1/Nf), hyper- 
geometric series give multiple zeta values (MZVs), such as ^(5, 3) = J2m>n>o ^/m^n^, 
in anomalous dimensions. This irreducible MZV already occurs at 0{1/Nf) in the e- 
expansion of the multiloop generator G{e,ne) in (|T^. 

Even the handling of the vertex renormalization (|16|), at large n, is nontrivial. For 
the n-loop term in the scalar correlator, one needs to expand the F functions and ra- 
tionals of to order and then make a further Laurent expansion of the in- 

tegrand. After integration, this three-fold series gets multiplied by J2k>o'^''^''~^j from 
the one-loop term. The fifth and final series is the most potent: one must multiply by 
{fi'^/Q^y = '}2k>Q{.LeY / k\ with the resulting complicated dependence on L = log(/i^/Q^) 
required to cancel nonlocal terms in the Laurent expansion of the multiloop diagrams. 
Thus the innocent-looking integral (p!6| ) generates a combinatoric plethora of products of 
rationals, zeta- values, powers of 1/e, powers of the coupling, and powers of logs. And thus 
far we speak only of the analytically tractable term, free of the MZVs that occur in true 
multiloop diagrams. 
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Next one sees, in Ib, the astounding interconnectedness of perturbative quantum field 
theory: there is a vast conspiracy, on a scale that would be ludicrous in human affairs, 
between analytically nontrivial n-loop integrals in ([T7|) , where G{e,ne) entails the ^3^2 
hypergeometric series of and the 5-fold series from the intricate combinatorical pro- 



cessing of g{e) by vertex renormalization. Everyone knows what this conspiracy must 
achieve: total cancellation of logs from the singular terms, as required by the locality of 
counterterms. How, one asks, is this conspiracy coordinated? Its success is not in doubt: 
field theory cannot fail. Yet nothing, prior to this work, appeared to indicate the analytical 
mechanism. 

Clearly the analysis of the vector case cannot solve this problem. Both of its key 
assumptions are now vitiated. First, the n-loop terms in (|T7p have a form that does 
not conform to (^); secondly, Laurent expansion of either ([161) or (|1^ generates a 



singularity at n loops, in defiance of the restriction < n in (|g). It was the Ward identity 
Zi = Z2 of QED, and hence of QCD at large Nf, that protected us from these eventualities 
in the vector case, by giving one less factor of 1/e, and by allowing an Ansatz (|^), with 
no singularity at n = 1. 



2.4 Reconciliation in the scalar case 

How is progress possible in this complex scalar case? By virtue of our recent finding that 

G{e,e)=g{e) (20) 

which expresses the remarkable fact that analytic continuation to n = 1, of the hypergeo- 
metric result for inserting chains ofn — 2 one-loop diagrams in a pair of two- loop skeletons, 
gives the anomalous dimension, to all orders in the coupling at large Nf. To prove (PU]), 
one must first show that the irreducibly hypergeometric terms in G{6, 6) vanish a.t 6 = e, 
hence cancelling the new singularity at n = 1 in (|T^), which was not encountered in the 
Ward-protected vector analysis, based on (H). Thanks to the systematic hypergeometric 



methods of ^^^^ is now possible. The surviving terms in the analytic continuation to 
n = 1 then give F functions, multiplied by a very complicated rational function of e and 
6. At 6 = e, the match of G{e,e) to the critical exponent g{e) is perfect. 

This enables us to organize both the combinatorics and the analysis, by writing 

G{e, S) = ^(e)^aM) {1 + {S-e) [G^(e) + Gn{S) + eSGs{e, S)]} (21) 
with a prefactor specified by 

Gi{e,6) ^ (^(/(^))'^')' (22) 
.^ _ F(l + ^) F(l-^ + e)F(l-2e)F(l-e) 
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and functions of e and 6 that we eventually obtain, by hypergeometric analysis, as 



GE{e) = e{l-2e) (24) 
^^^^^ - 1^ " 2^ + (P - (1 - 5)2)2 (25) 

= E - '^-V-' - I E C(2/ + 1)/(1 - A-'W'-' (26) 

fc>o ;>o 

with (p^) showing the renormalon structure and (p6|) giving the Taylor expansion about 



5 = 0, thanks to (|TT],|T2D. As in the vector case [|T4|, [T8| no even zeta value can occur in the 
expansion of the renormalon contributions, since = I{—S) is conformally invariant. 
In the new scalar result (p6|) one has the further simplification that odd zeta values occur 
only in odd Taylor coefficients. The vastly more complicated residuum, G3, which we also 
determined completely, makes no contribution to either the anomalous dimension or the 
finite part, since it is multiplied by a factor that vanishes at 6 = e, 6 = 0, and e = 0, 
corresponding to n = 1, n = 0, and n — 00. After extracting unity, Ge, and Gd, one 
may simply throw away whatever remains in G3. We did this the hard way, by calculating 
everything, exactly. 

The highly coordinated conspiracy (pO|) is signalled by the leading unity in the braces 
of (^11) , and the triviality of Gi{e,e) = 1. Moreover, the remarkable combination of F 
functions in (p3|) gives G2{s,e) = G2{s,0) = 6*2(0,5) = 1, which means that we may 
set G2 to unity, without the shghtest effect on the outcome, for the same reason that we 
may set G3 to zero. Next we note that (p^ ) precisely cancels the rational denominator 
of the one- loop term (0): this renormalon- free term is as simple as one could ever hope. 
Inspecting the analytic structure of the nontrivial renormalon contribution (^) we see 
single poles at 5 = 1 and 6 = 2. The pole at 5 = 2 corresponds to the appearance 
of a gluon-condensate contribution of order m?'Q'^{G^^yG^'^) /Q^ in Ilg. This was to be 
expected, by analogy with the vector case. 

The totally new feature is the pole at 5 = 1, reflecting the ambiguity of the constant 
term in the correlator, which is inflnite, though formally proportional to m"^, in pertur- 
bation theory, while current algebra relates it to {m^xjj). It is wonderful that an analysis 
of finite parts of massless diagrams at large loop numbers can so powerfully remind one 
of what one knows from massive diagrams 0] at low loop numbers, namely that the UV 
physics of the second subtraction in the dispersion relation is as profound as that of the 
first. It makes no sense to say that at large we can forget about the second UV sub- 
traction, because m'^/Q'^ is small. If we make believe that we can, the 5 = 1 infrared 
(IR) renormalon of massless diagrams will remind us of the ultraviolet physics of massive 
diagrams. It was, of course, the Ward identity that protected us from this consideration 
in the vector case. 



Proceeding, we see that the closed forms (p^ , |25| ) lead to renormalized contributions 
that are speedily found, to very high orders, by the combinatorics (||,|^ that served in the 
vector analysis of |]T3|, where results were obtained up to 20 loops, analytically, and up to 



100 loops, numerically. It remains, however, to resum the leading term in (pi]). We may 
set G2 = 1, and write l/Gi = (1 — Gi)/Gi + 1, with only the final unity requiring further 
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attention. That is now feasible, since the formal transformation 

gives a rather simple Laurent series, when multiplied by g{e) = J2n>o9n^"' ^- Finally, we 
expand l/{e — x) in x/e and formally integrate (p^. 

Tidying up, we see that Ia + Ib is equivalent to J a + Jb, as far as pole terms and finite 
parts are concerned, where 

E^P^-E^^ (28) 

combines I a with the recalcitrant leading term from Ib, and the remainder of Ib is equiv- 
alent to 

j«(M)..s(.)i:f^)""^^- (29) 

^i\b-ej Gi{e,ne)n 

where 

Lsis, 5) = + £(1 _ 2s) + GoiS) (30) 

e — 



is easily processed by (|g,0), using ([19|j2^j25|) . It is important to note that J a is no longer 
minimal: the first sum combines pure pole terms of I a, with weight l/(n — 1), and the pole 
terms from multiplication of g{e) by (pT]), with weight — 1/n. Then a regular part appears, 
in the second series of J a- The overall 1/e singularity, from the one-loop diagram, means 
that the k = term in (|28D becomes singular, and the k = 1 term becomes finite. This 
will be seen to be crucial, in the following analysis. 



2.5 Critical behaviour of the scalar correlator 

Now we analyze the so-called scalar-scalar anomalous dimension, jss, defined by |T|] 

^ ^ + /5(«s)ttt^ + Us = [m(/i2)]2g2 + o(^]] (31) 



log /i^ a log tts / y \Q / J 

The presence of an 0{m?' /Q"^) term in the braces on the r.h.s. of the renormalization-group 
equation ( PT| ) reminds us that two subtractions |^ are required for the scalar correlator. 



Here we are concerned with the first, which produces a scale dependence described by 7^5. 
If one supposes that the second might be forgotten at large Q^, the unity of field theory 
soon corrects one: we have already seen that the finite parts of massless diagrams at large 
loop numbers are profoundly aware of the Oijn^ /Q"^) UV physics in (|3l[) , via the 5 = 1 IR 
renormalon in (|25D, at m^/Q"^ = 0. 

Working to leading order in l/Q"^, and next-to-leading order in 1/Nf, we write 

= [mif,')]'Q'dF(^-2L-A + ^^H{L,b) + 0il/Q') + 0il/N])^ (32) 

iss = rf^(-2 + ^M6)+0(l/iVj)^ (33) 
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where b = TpHfas/Sir is the leading term of the 4- dimensional beta function at large Nf, 
i.e. the value of e such that the critical dimension is 4 — 2e. Hence the renormalization- 
group equation (0) simplifies to 

+ ^'^) ^H{L, e) + 4(L + 2)eg{e) = eh{e) (34) 

where L = log(/i^/Q^). Our aim is to determine h{e) to all orders. Then the dependence 
of H{L, b) on L is completely determined by H{0, b). 



We stress the underlying principle of this work by using the argument e in (|3^) , where 
a reader more used to 4-dimensional perturbation theory might reasonably expect to see 
us choose b = TpNfas/3Tc. The reason should be clear: the large- A^/ beta function is b — e 
in (i = 4 — 25 dimensions, where e need not be small. By working near the critical point, 
b = e, without assuming that d is near 4, we bypass perturbation theory. 

We have remarked that true anomalous dimensions differ from (^) only by rational 
functions of d, at order l/Nj. Of course the beta function is not such an object, since 
it vanishes at the fixed point, by definition. The reason why the 0{1/Nf) corrections to 
the 4-dimensional QED beta function are given by the integral of (|]) is clear: the physics 
resides in the critical exponent that is the derivative of this integral. In the scalar case, 
the critical exponent is the derivative of eh{e) + 4:g{e). With some ingenuity, one may 
obtain it by careful parsing of (|28|J30|), which yield 



1 d f eh{e) + Ag{€)\ (rf - 3) 



2 



g{e) de\ 4 / 
Here, in short order, is the proof of this fine parabola: 



-2. (35) 



1 {f{e)-l)/e + e{l-2e) + 2 \ _ l-f{e)-e' _ {d - 3)' ^ ^^^^ 



d-3[ fie) J efis) 

with a seemingly troublesome singularity at d = 3 turning out to give a harmless minimum. 
The prefactor l/{d — 3) = 1/(1 — 2e) comes from the one-loop result ([T5|) . The first term in 



braces comes from the first series in Ja; the remainder from setting 6 = in Lb- The weight 
of the series is doubled by the shift 1/6{1 — 2e) = l/£: + 2/(l — 2^); the fourth term follows 



from using ([T0|) in (pSj), which gives GniO) = 2. By taking the derivative of eh{e) +4(7(6:), 
we remove the nonminimal series in J^, whose weight, 1/n, is incommensurate with (|]), 
which we use for the Jb terms. The rational function j[e) enters via (^2]). Taking the 
exact one- loop vector result from (|13|), we obtain parabola (|35D . 

Recall what makes this possible: the circumstance (|20|) that analytic continuation 
of the hypergeometric series of [0 to n = 1, i.e. to minus one insertions, reproduces 



the large- iVj anomalous mass dimension to all orders in the coupling. We began with 
two bad problems, for which the vector analysis gave no preparation: first we had an 
extra singularity, vitiating the combinatorics; secondly we had incommensurate weights 
for mass renormalization and those terms that we could handle. The beauty of ( [2(]| ) is that 
it enabled us to solve both problems at the same time: transferring the combinatorically 
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recalcitrant term to join the mass renormalization in J a, we obtained the desired weight 
as l/(n — 1) — 1/n = — 1), in conformity with (^. There was a price to pay: 

this transfer took a nonminimal term with it. But that was no problem: we knew from 
critical phenomena that we must take a derivative, so as to obtain a physically significant 
exponent. In that derivative we must include the precise multiple of g{e) that kills the 
nonminimal terms. 

The road from analysis of the Saalschutzian ^3^2 series of multiloop diagrams, to the 
simple parabola (|35D , was a long one. Along it there was a narrow bridge: G{e,e) = g{e). 
We have found no other route. 



2.6 Analytical results at large Nj 

Recalling that h = TpNfas/3iT, we obtain 

755(«s) = dp (^-2 + ^ E {^^^f^y ' + 0{l/Nf)}^ (37) 

where h{e) = J2n>i hnB"'~'^. Working to merely 4 loops, we immediately find 

1 e; 162^ 
/i2 = -y, h = 9, /,4 = -18C(3) + ^, (38) 

in agreement with [|l|. The development 

h = -27C(4) + ^C(3) + ^ 

hr = -90C(6) + 30C(5) + ^C(4)-18[C(3)p + ^C(3) + ^ 
h = -162C(7) + ^C(6) + ^C(5) - 54C(4)C(3) + ^C(4) 

was obtained , using 



hn+l + 4^n+l _ 4gf„,2 - 45'n.-,i - 

4 ~ 3n 

which solves (|35|) , with 



n>l 



2" 2 



(39) 



exp|E '' '^V (O^M (40) 



obtained from (0). By this means analytical results to 20 loops, and numerical results 
to 100 loops, are readily obtainable. Since g{e) is finite for e < 5/2, so is h{e). Hence 
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the coefficients decrease rather rapidly, with hn = 0((2/5)"). For example, we found that 
/lie ~ 1-867 X 10-^ 

This convergence is in marked contrast with the behaviour of the finite parts. Writing 



H{L,b) = J2n>i Hn{L)b'^ ^, we used (p8| , |30D to derive the all-orders solution 



nin-l)Hr,iL) = - 4(L + 2^) + - 9(-irZ}„(L) (41) 

5^D„(L)57d = {l + 5GD(5)}exp((L + 5/3)5) (42) 

n 

of the renormalization-group equation (^4]). Two good checks are provided by the vanishing 
of the r.h.s. of (|T|) at n = and n = 1, using gi = 9/4, /i2 = 4^(2 = -15/2 and Gd(0) = 2. 
Two stronger checks are provided by known results at two and three |jl| loops. We 
obtain from ( ^T]) the entirety of the former and the large-A^/ terms of the latter, for all 
/iVQ2. At At^ = (5^ we write if„ = H„{0). In [g, we find 



H. = 9 



(43) 



at /i^ = g2_ 

In order to extend explicit results to high loop numbers, and to analyze their UV and 
IR renormalon content, we separate Gd(5) = G^{6) + G^{6) into 

and expand in S. (Note however that the Taylor expansion of the total renormalon con- 
tribution is already given with great economy by (^6]), with odd zeta values in odd Taylor 
coefficients.) Computation of analytical results to 20 loops takes seconds, for any value of 
L = log(/i^/Q^). The analytical development to merely 8 loops is given by 

X 27^, , 157^, , 499069 
= 90C(5) - -C(4) + — C(3) - 

- -^C(5)-^C(4)-^C(3) + ^^ 
5 5 ^ 4 ^ 60 ^ 1728 

= 1701C(7)-15C(6) + ^C(5) + ^C(4)-31C(3)l^ + 51^C(3)-?i^5£15? 

b sv ; ^ ^ 10 ^ ^ 40 ^ Lsv j\ -r sv ; 155520 

Hr - -^C(7) - i^C(6) - ^C(5) - ^C(4)C(3) + f C(4) 

_ ^\C(3)? - ^^C(3) + ^^^^5753995 
)\ 189 ^ 40824 

, 1323,, , 3967083, 3635,, , 27,, , 7017943,, , 
Hs = 68850C(9) - ^C(8) + ^^C(7) + ^C(6) - yC(5)C(3) + —^d^) 

639,,,,,, , 12973,, , 727,,, ,,2 28819423,, , 315995418895 

C 4 C 3 + C 4 + C 3 f + C 3 

ggsv AV ggg sv ^^112^^^ ^ 72576 ^ 1492992 
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2.7 The asymptotics of naive nonabelianization 



We reiterate the often rehearsed (yet never properly justified) argument for Borel resumma- 
bihty of large- A^j singularities at 5 < 0, in the case of QCD. At truly large Nf, they give 
sign-constant series. Now we imagine that, in some vagueQ sense, the gluon loops "follow" 
the quark loops. Hence we perform the naive nonabelianization []T5| Nf — >■ Nf — ^^^^ 
(iVc = 3), which gives b -poa,{fj,'^)/7i, with po = {IIN^ - 2Nf)/12 giving the relative 
contributions of gluons and quarks in the one-loop beta function of QCD. In the real 
word, with Nf < 6 active quarks, we have Po > and b = /3oas(/i)/7r ^ 1/ log((5^/A^), if 
we suppress large logarithms by renormalizing at fi^ = Q^. By this sleight of hand, the 
singularities at 5 < now give a sign-alternating asymptotic series that is resummable by 
Laplace transformation. A singularity at 5 > leads to an intrinsic ambiguity of order 
exp(— 5/6) = 0{Q^'^^) in a dimensionless correlator. In the vector case the dominant 
ambiguity, at 5 = 2, is associated with the gluon condensate. Here, in the scalar case 
with Us = 0{Q^), a singularity at 5 > leads to an ambiguity of order Q"^'"^^ . Hence the 
renormalon at 5 = 1 is in accord with the fact that the constant term in H5 is inaccessible 
to perturbation theory. Interestingly, our analysis does not distinguish the scalar from the 
pseudoscalar channel. In the latter case, current algebra relates the constant term in the 
correlator to the quark condensate (jniptl)). 

At any given order n > 1, we may separate if„ into 6 parts: 

= + + Ht^ + + + (46) 

where 

^Ms ^ n{hn+i - 4(L + 2)gn) + Ag^+i 
" n{n - 1) 

is the highly convergent part from MS renormalization and H^^^ comes from the lead- 



ing unity in the braces of (|42|) . The next term comes from (0), with resummable UV 
renormalon singularities at 5 < 0; final three terms come from the IR renormalons in (^5]). 

Table 1 shows the numerics of this breakdown, again at //^ = . We comment on 
each contribution in turn. 

1. The MS-specific contribution = 0((2/5)") is negligible at large n. This is 



because the critical exponent (19) is finite at all dimensions d> —1. The numerator 



of (^) shows that the same applies to the critical exponent in the vector case. 

2. The modest growth of H^^ comes from the choice iJ? = Q^, made so as to com- 
pare (H) with |ll. At /i^ = Q'^exp{-5/3), corresponding to QED MOM- or V- 
schemes subtraction of quark-loop insertions, these terms fall off like 1/n^. 

3. The series coming from H^~^ is now regarded as infrared- safe: one may resum it by 
Laplace transformation of the Borel transform (PD, obtaining 



n>l ^0 ^ 



^It must be vague. Unlike quark loops, gluon loops are not gauge invariant. 
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where b ^ 1/ log(/i^/A^) replaces —b, by naive nonabelianization. 

4. The total is dominated by the factorial growth of H^^ = 0{{n — 2)!). Indeed 
this dominant (5 = 1 IR renormalon even gives a reasonable account at very low 
orders. For example, it gives a fraction 73.0/84.5 = 86% of the 3-loop large- A^/ 
result of 1^. Clearly perturbation theory would be in bad shape, if this renormalon 
entered phenomenology. Fortunately it does not. Rather it is a reminder, from the 
perturbative sector, of the long-distance physics in 115(0), which current algebra 
relates to quark condensates. In the sum-rule analysis of 0, it was nullified by 
taking a twice-subtracted dispersion relation. In Higgs decay and spectral function 
of the QCD sum rules, it is nullified by taking the imaginary part, in the physical 
region s = —Q^ > 0. 

As previously remarked, one must set aside the temptation to divide Us by and 
then differentiate w.r.t. Q^, so as to remove 755 from (|3T|). That would not remove 
the 6 = 1 renormalon, since it would leave the kinematic singularity Ils{0)/Q'^. Then 
long-distance physics would make the asymptotic perturbation series expire sooner 
than needs be, since it is always faithful to the motto dulce et decorum est pro 
Wilson mori: it's OK to explode in accordance with Wilson's OPE. Indeed, it often 
signals impending doom at the 2-loop level . 



5. More long-distance physics resides in H^^ = 0{{n — 2)!/2"), which signals the pres- 
ence of the gluon condensate in Wilson's scheme of things, as articulated in QCD 
by p6|. This 5 = 2 IR renormalon is also absent from the imaginary part, to lead- 



ing order in 1/s. When taking an imaginary part, in the physical region, one kills 
any single-pole renormalon, and turns a double-pole into a single pole, since the 
imaginary part of ( ^2]) acquires a factor sin(7r(5). The quark and gluon dimension-4 
operators appear in the energy momentum tensor and are hence renormalization- 
group invariants, to leading order. Thus they are absent from the imaginary part, 
at order l/f3o, in the high-energy limit. Since there is no matrix element to absorb 
renormalons at 5 = 1 or 5 = 2, these, and only these, appear as single poles in (^Sf). 



6. The only IR (i.e. unresummable) large-/3o renormalons that appear in the imaginary 
part at large s, are those in at 5 > 2. These correspond to long-distance 

physics in matrix elements of operators Ok with dimensions dk > Q in the OPE. The 
resultant ambiguity in Im II^ is of order m?{Ok) / s'^*'~^ . 



2.8 Analysis of the imaginary part at large Nj 

Table 1 might appear alarming. A far happier picture emerges in Table 2, where we 
analyze a physical quantity, namely the imaginary part Im II^ = 2ttsRs{1 + 0(l/s)) at 
—Q^ = s = w"^. For w = Mh, this contains the radiative corrections to decay of a Higgs 
boson of mass Mfj into a quark- ant iquark pair, ignoring terms of order {m{M'jj) / MhY . 
Now we are dealing with a multiplicatively renormalized quantity: the explicit dependence 
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of Rs on is cancelled by its implicit dependence, via Q;s(/i^) and m{ji^), giving 



+ /5(«s)tti + 27m(as) 1 i?5 = (49) 



\(91og/i^ SlogOs 

with appearing because Rs contains two powers of the renormalized mass in the Born 
approximation. We should, however, set iJ? = 0{s), to suppress large logarithms. Here, 
we set fi"^ = s = and obtain 



Rs = 3[m{w)r f 1 + - E bs-' {Sn + 0{1/Nf)}] 

\ ^ n>l J 



(50) 



with bs = —Nfasiw"^) /6tt, which is replaced by /3oa<i(w^)/vr ^ l/log(w^/A^) in naive 
nonabelianization. We have 5*2 = 17/3, with the large- A^/ term giving the total radiative 
correction at 2 loops. At n > 2 loops, we obtain the leading term at large- A^j as 

Sn = A„ + A„ (51) 

A,^?tDli. ,52) 

n-1 gi 

= 2eM55/3) {I + 6G45) + 6GM}"^^- 2 (53) 

where Qn is the n-loop term in the large-A^j result (^D|) for the anomalous quark-mass 
dimension. The separation Sn = An + A„ into anomalous- dimension and renormalon 
contributions will be used in Sec. 3, where we shall retain only the latter in large-/?o 
approximations. Here, we retain both terms, so as to present the exact analytical results 
at large- A^j. 

We note that the high-energy imaginary part Rs receives powers of vr^ from two sources: 
from the even zeta values {({2k) \ k > 1} of euclidean analysis, and also from the analytic 
continuation of logarithms to the physical region. At large- A^y, the separation of these 
two effects is particularly clean: even zeta values occur only in the anomalous-dimension 
contributions An, powers of vr^ from analytic continuation result only from the factor 
sm{'ir6) /tcS in the Borel transform ( ^3]) of the renormalon contributions A„. We write the 
An terms in braces in the following large- A^/ results. Up to 4 loops we obtain 

f 51 22 17 

Ss - m-4C(3)-l.^ + ^^-4C(3)-l.^ + ^ (55) 
l36j ^ 3 18 ^^^34 ^ ' 

r4^, , 83 1 40^, , 22 2 3940 , 22 2 15511 , , 

^4 = -C 3 > C 3 TT^ + = -12c 3 TT^ + 56 

13^^^ 108] 3^^^ 9 81 ' 9 324 ^ ' 

with totals in agreement the with the large- A^j 3-loop 0] and 4-loop |l| terms. At 5 and 
6 loops our new results 

3.,.. 5._. 65- 
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MC(5) + 4C(3).^-Mc(3) + i.'-^.^ + 5|II (57) 

12.... .... 7.,„, 451- 



= ^. + A. = {_C(5)-C(4)-^C(3)-^} 



sv ; ^ 3 sv ; 27 ^ 15 81 243 ^ ^ 

were found by expanding (|40|) , to obtain A„ in the braces, and the new all-orders renor- 
malon results ( P| , ^5D , to obtain A„ via (|^). This method may easily be continued up to 
20 loops, analytically, and up to 100 loops, numerically. 

The corresponding vector quantities in electron-positron annihilation, at large Nj, 
come from the old result (|^) of |T^, which immediately gives 



^-^K ^ 8exp(5V3) /_ (-1)^ ^ (-1)^ ] sin(7r^) 

tAn-2y. 3(i-5)(2-5) 1^ toik + sy^t2ik-^r) ^5 ^ ^ 

At 6 loops the vector coefficient has grown by an order of magnitude and changed sign, 
with Vq/V2 ~ —11. By contrast, 8^/82 ~ 0.5 in the scalar channel. This remarkable 
postponement of factorial growth, at large- A^/, does not depend on a cancellation between 
anomalous-dimension and renormalon terms; rather it reflects cancellations between the 
renormalons themselves, with A6/A2 ~ 0.4 showing no sign of the growth that had become 
clear at 6 loops in the vector channel. In this respect, the (pseudo-) scalar channel is better 



behaved than the (axial-)vector channel, at high energy and large Nf^ despite warnings |52] 
that might suggest an opposite situation. 



2.9 Postponed factorial growth 

Considerable interest attaches to the numerics of Table 2, where it will be seen than S'„ 
is amazingly well-behaved for n < 7, while Hn in Table 1 had already gone haywire at 
n = 3. We have discovered a plateau of tranquility at loop numbers n = 2, 3, 4, 5, 6, in the 
large- iV/ terms of the imaginary part of the scalar correlator. 

The behaviour at n > 7 in Table 2 is fairly clear. 

1. The MS term 8^ falls of like (2/5)". 

2. The leading unity of the braces of (^) gives power growth of 8^\ modulated by a 
sine. 

3. Eventually, the UV renormalons at 5 < take over, giving an alternating series for 
8^\ now that we have naively nonabelianized. But they takes ages to get going: 
even at n = 8 they have not yet overtaken the humble 6 = term. 

4. The term 5*^^^ from 5 = 1 is no longer a renormalon; the single pole has been cancelled 
by sin(7r5). 

5. The same applies to the 6 = 2 term 8^\ except that it is, on average, smaller than 
the 6=1 term, by a factor of order 1/2", after allowing for the sinusoidal oscillations. 
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6. The unresummable series with the coefficients from the IR renormalons at 5 > 3 
grows factorially, eventually. However it is suppressed by a factor of order 1/3", in 
comparison with the UV renormalons. 

The staggering feature is the tally, Sn- For n > 7, it behaves like the wild animal 
that it truly is; for n < 7, all is sweetness and light. Presented with only the analytical 
expressions for S'„ at n = 2, 3, 4, 5, 6, one would not have the slightest inkling of what is in 
store. Conversely, one may say that the large- iVj renormalons show themselves mercifully 
late, in this physical quantity. 

Note that these statements are entirely dependent on the scheme chosen (for a dis- 
cussion of the scheme- dependence of renormalon contributions see e.g. Refs. |2^). As 
will be clear from the discussions of Sec. 3.5 and Table 5, the V^-scheme results are in 
fact a much better indicator of the eventual asymptotics of the coefficients. Nonetheless 
given the widespread use of the MS scheme in the literature we initially focus here on the 
asymptotic behaviour in that scheme. 

Now consider the case where s = w"^ with w ^ 2 GeV, as in strange-quark mass 
extraction With Pq = 9/4, the expansion parameter Poasiw)/^ ~ l/log(i(7^//i^) ~ 



0.2 is now uncomfortably large. Hence the tranquil plateau at n < 7 loops is good news 
for the later analysis in 0, where 4- loop perturbative QCD was used, on both sides of the 
sum rule. The result was significantly different from an earlier 3-loop analysis in 0], by 
the same group, using a different truncation procedure (compare with the similar 3-loop 
studies of Ref.lSl). 



The villains, which might have been waiting just round the corner, were the renor- 
malons. Might a large 5-loop term significantly change the 4-loop result? 

To date, we know of only one analytical technique for estimating such effects on the 
basis of genuinely new calculation, instead of reshuffiing old input: naive nonabelianization 
of the large- A^j terms, which we have computed in the demanding scalar channel, at some 
cost of labour. The good news is that we found nothing alarming at n = 5 loops. The 
even better news is that all seems well at n = 6 loops. Only at n = 7 does the inevitable 
growth show signs of commencing. Hence the best indicator that we can compute suggests 
that the perturbative part of the strange-quark-mass extraction in is in fine shape. 

Indeed, the results of Table 2 indicate that the contribution of the 5-loop coefficient 
might be not crucial and that the corresponding perturbative series should be truncated 
at the 6-loop level in accordance with the common practice in treating the predictions 
of asymptotic perturbative expansions, which presumes their truncation at the minimal 
term. 

This fits nicely with the claim in that the condensate contributions are also under 
control. Had these contributions been substantial, our discovery of perturbative tranquility 
at large Nj would have been rather puzzling; now it may be taken as gratifying evidence 
of the depth to which the OPE connects ultraviolet and infrared physics. 
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2.10 Euclidean analysis at large Nf 



In the vector channel we may express the results of high-energy perturbation theory in 
two ways: either in terms of R{s), in the physical region, or in terms the Adler function 

in the euclidean region. Here, R{s) = 1 + a^/vr + O(a^) gives the high-energy radiative 
corrections to the parton model in electron-positron annihilation and hence -D(Q^) gives 
the corresponding radiative corrections to the derivative of the polarization function of 
the vector correlator. 

In both cases, we ignore quark masses, so the transformation between one set of radia- 
tive corrections and the other is generated, to all orders, by the following integral 



with the expansion of 



TlS 



sin(7r5) 



1 + ^ (2 - 4^-^) ({2k)S^'' (62) 



fc>0 



telling one precisely how to remove from the imaginary part all and only those powers of 
TT^ that came from analytic continuation of logarithms. 

In the scalar channel, we are confronted by a choice of euclideanizations of the high- 
energy imaginary part. The most prudent choice would appear to be the dispersion relation 



for the seconc? derivative of the scalar correlator [|2|, which is multiplicatively renormalized 
for all values of m'^/Q'^ and is hence free of the IR renormalon at 5 = 1 in (^). Ignoring 



terms of order m'^/Q^, this amounts to the euclideanization 

Ds{Q) = 2Q (s + Q^y^ ^ I t — n — + ^*^"^'') 



of the radiative corrections to the imaginary part at high energy. However, a merely 
mathematical analogy with the vector case might lead one to consider the construct 

corresponding to a dispersion relation for the first derivative of Ils{Q'^)/Q'^ ■ Here one 
expects the asymptotic perturbation series to destroy itself earlier, leaving an ambiguity 
of order A^/Q^ that reflects the failure to remove the infinities in H5'(0). Therefore, this 
ambiguity has a perturbative origin. 

We shall show that at large- A^/ the perturbative series for the nonstandard euclideaniza- 
tion Ds is indeed worse behaved than the twice-differentiated euclideanization Ds, in ac- 
cord with the expectation from the OPE. It might therefore be expected that we shall 
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proceed in Sec. 3 only with the safer alternative Ds, as used in QCD sum-rules 02 



In fact, we shall need both constructs, so as to study the logarithmic derivative of Ds , 
which is the euclidean analog of the considered in the case of Higgs decay quantity (see 



Refs.lBI, O, |56ll). In our case it can be defined as 



which satisfies a renormalization-group equation 

2 + ]Rd = (66) 

y o log fi^ o log as J 

that is free of the anomalous quark- mass dimension and hence suitable for analysis by the 
method of effective charges (see also Ref.||58||), scheme-invariant perturbation theory 



59| , commensurate scale relations and the standard PMS approach |l6T| 



Here we assemble everything that is known about the MS perturbation series of the 
euclidean constructs ( |63| , |64| , |65D . It is convenient to begin with 

DsiQ') = 3[m{Q')r(l + Y.dn(^^Y] (67) 



n>0 



TT 



10801 39^,^, /65 2^,^,\^^ 



V 486 9 

obtained by removing terms involving vr^ in the expansion 

'as{w) 



Rs{w') = 3[m{w')r fl + E^n f- 

\ n>0 \ 



TC 



(71) 



of the imaginary part, given to 4 loops in [|T|. To effect the inverse transformation, to 5 
loops, one may use the fixed-order perturbative expansion in the minkowskian region 

si = rfi (72) 

52 = ^2 - 7o(/5o + 27o)7rV3 (73) 

53 = 4- [rfi(/?o + 7o)(/5o + 27o)+A7o + 27i(/3o + 27o)]7rV3 (74) 

54 = ^4- [rf2(/?o + 7o)(3/9o + 27o) +c/i/?i(5/?o + 67o)/2 + 4di7i(/3o + 7o) 

+ /527o + 271 (A + 7i) + 72(3/5o + 47o)]7rV3 

+ 7o(/?o + 7o)(/5o + 27o)(3/?o + 27o)vrV30 (75) 



19 



where the relation between S4 and was derived in Ref.pS] and 

n+l 



dlogm 



^log/^' n>0 



(76) 



gives the expansion of the anomalous quark-mass dimension, in the same manner that (|I|) 
gives the expansion of the beta function for the scale dependence of the coupling. Both 
expansions are known to 4 loops. The coefficients of 7.^ are |]11], O]: 



7o 
71 

72 

73 



1 

1 

16 
1 

64 
1 

256 



202 



3 

1249 
4603055 



20,, ■ 
—Nf 
9 ^ 



/2216 160 



V 



140 



162 27 
/91723 34192 



+ C(3) ]Nf 

27 3 V 81 ■ 

135680,, , 

C(3) - 8800C(5) 



+ 



+ 



V 27 
/5242 



+ 



9 



C(3) - 880C(4) - 



18400 



V 243 



+ 



^00 
"9" 



C(5) 



9 

^332 
243 



64 
27 



C(3) iV: 
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while those of beta function are ||13|| : 

2 

3" 

102 - 
■2857 



1 

4 
1 

16 
1 



64 
1 







38 




y 


Nf 


5033 



-N 



f 



256 



6 



325 ,,2 
54 L 

1078361 6508,,^,,,, 
+ -,^a3)]Nf 



162 
1093 



27 



/50065 6472,, A,. 

+ + C 3 ]N^f + 

V 162 81 V ^ 729 



which will likewise be needed in our analysis. 
The results at large- A^/ are 

dn-i = (-iV//6)"-' 

n>2(^-2)! 



A^ + A,, + Oil/Nf) 
2exp(55/3) {1 + + 6G+{6)} - 2 



with the Borel transform (|86D giving 

A5 : 



, 100,, , 64877 
-60C(5)--C(3) + — 



~ . ^ 2000,, , 244871 

Ae = -400C(5) - — -C(3) + 



(77) 
(78) 

(79) 



(80) 

(81) 
(82) 
(83) 

(84) 



(85) 
(86) 



^7) 



27 



243 
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at 5 and 6 loops. 

We expect the factorial growth of (pSj), in Ds, to be more drastic than that of 

Sn-i = (-iVy/6)"-2 {An + A„ + 0{1/Nf)) (89) 

in the imaginary part Rs, since the latter is free of the spurious 6=1 renormalon that 
afflicts the former. Table 3 emphatically confirms this expectation at large Nf, where one 
sees that A6/A2 ~ 69 is two orders of magnitude larger than A6/A2 ~ 0.40. 

In general, one expects that the asymptotic structure of perturbation theory expansions 
will differ for the physical quantity Rs{s) and the euchdeanization Ds{Q^)- To construct 
the latter from the former one takes the renormalizat ion-group determined powers of the 
minkowski logarithm log(/i^/t(7^) and performs the transformation 

- 4^-'=) C(2A;) (90) 

on even powers. If the imaginary part is fairly well behaved, as at large Nf, it is unlikely 
that its euchdeanization will be so. Indeed, the factorial growth of the r.h.s. of Eq.(|90D 
should restore the factorial growth of the perturbative series for the euclidean quantity Ds, 
expected from general grounds of quantum field theory. The large- A^/ analysis of Table 3 
suggests that the imaginary part Rs is rather well behaved for n < 7 loops, with the far 
worse behaviour of Ds (and thus A„) resulting from its renormalon at 5 = 1, which was 
suppressed by the sine function in (^). 

Next we consider the more prudent euchdeanization (§^), with a perturbation series 
Ds{Q') = 3[m{Q')r(l + Y.^dn(^^^J'^ (91) 

n-l 

dn = dn - 2-fn-l - "^{kPn-k-l + '2-fn-k-l)dk (92) 

k=l 

The large- iVj results for Ds are obtained from 

dn-i = (-iV^/6)"-2(A, + A,, + 0(l/iVy)) (93) 
A, = A„ - (n - 2)A„^i (94) 



at n > 2 loops, with the 6 = 1 renormalon removed by the combination (Q). At 2 loops, 
we have A2 = A„ — 2 = 16/3. Table 3 shows that the factorial growth in Ds at large Nf 
is milder than in Ds, but more severe than in the imaginary part Rs- 



From (^^2D we obtain the first 5 terms in the expansion 



Rd{Q') = E-nf^^V"^ (95) 



TT 



n>0 

ro = 70 = 1 (96) 
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n 


= 71 + i^l3odi 




(97) 




= 72 + (^od2 - ]^f^odl + ^/?ic?i 




(98) 




3 13 

= 73 + 2/^o4 + 2-^0^?? - -^Podid2 


+ M2 - \mi + \Mi 


(99) 




= 74 + 2/?o4 + 2podld2 - 2Podids - ^Podt - Ml 






3 13 

+ + -Ml - ^Pidid2 - 


]^Ml + M2 + ]^Mi 


(100) 



where the 5-loop coefficient d^ is unknown, while the estimates for the 5-loop anomalous 
quark mass dimension are known from the results of application of the Fade resummation 
method ( note, however, that the analytical calculation of order 0(l/A^|)-corrections 
to 7m(as) [0 indicate that the latter ones should be refined). At large Nf we simply 
obtain 

r„ = {-Nf/QT (^A„+i + Oil/Nj)) (101) 

for > 0, with no contribution from the anomalous quark-mass dimension beyond the 
one-loop result tq = 70 = 1- It follow that at large-A^/ the perturbation series for i?^) 
explodes as violently as that for Dg- 



3 Subtleties of the naive nonabelianization 



In the previous section we mainly concentrated on the analysis of large Nf perturba- 
tive results, for the different theoretical quantities, related to the correlator of quark 
scalar currents. However, as was already explained in Sec. 2. 7, it is of definite interest to 
study the truncated and Borel-resummed perturbative series within the framework of the 
NNA Ansatz, which is postulated by applying the substitution Nf ^ Nf — llNc/2=—6Po 
(where Pq is the ffist coefficient of the QCD /^-function, defined by Eq.(l)), supplemented 
by retaining the leading terms in powers of f3o in the reorganized perturbative series. This 
procedure enables one to transform large- A^/ results, which are related to QED (note that 
in QED is proportional to Nf), to the nonabelian case of QCD. 

In this section we shall study a number of theoretical issues related to the application 
of the NNA approximation in the scalar channel. In particular, we shall concentrate on 
obtaining estimates of uncalculated higher-order terms in perturbative series for quantities, 
related to the correlator of quark scalar currents both in the euclidean and minkowskian 
regions. We shall also formulate different procedures for the resummation of the large 
minkowskian vr^-terms, within the framework of the NNA approach. 



3.1 Estimates of the higher order corrections in the MS-scheme. 

We begin by considering the expressions 

^NNA ^ /3-iA„+i (102) 
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^NNA ^ ^-iA„+i (103) 



for the coefficients (i„, s„ and (i„ in Ds-, Rs and at n + 1 loops in the MS-scheme. 
These are obtained by the naive nonabehanization of the terms in (^,^,^) that do not 
involve the large- A^/ anomalous quark- mass dimension contribution An, which is known 
to be unamenable to naive nonabehanization in general (for detailed discussions, related 
to the deep- inelastic scattering anomalous dimensions see the works of Ref . [|6^ ) , and in 
any case is small at 5 loops, since it falls off like (2/5)". 

Table 4 shows that all three NNA estimators give the correct sign and order of magni- 
tude at 3 and 4 loops. At 4 loops, the success of d-^ is rather remarkable, since we are 
using only the {Nf — 33/2)^ approximation to a quadratic to estimate the full result. In 
particular, at Nf = 5, as in Higgs decay, differs from the NNA estimate by only 8%. 
It is significant that this success of NNA occurs in the safer euclideanization Ds, which 
includes neither the tt^ terms of Rs nor the spurious renormalon of Ds- Accordingly we 
take 

3r = ,3„'(^ + fc(3)-60C(5)) (105) 

as our favoured NNA estimator. 

Note, that even if we choose an overall factor of two as the conservative uncertainty in 
the estimating power of the NNA procedure (which is motivated by inspecting the related 
numbers of Table 4 for Nf = 3), we arrive to the conclusion that the 5- loop perturbative 
approximation for Ds is really well-behaved. Indeed, taking Nf = 3 and as ~ 0.3 we 
obtain the following series 

Ds = l + 3.67(^^) +14.17(^^)' + 77.36(^^)' + 2x 1.26(^^y (106) 
= 1 + 0.350 + 0.129 + 0.067 + 0.0002 

with a rather small 5-loop term. Thus, the 4-loop extractions of the running mass rris 
of Ref.0, which is based on the consideration of the D^-function, indeed contains rather 
small perturbative QCD uncertainties due to the truncation of the corresponding series at 
the 4-loop level. 

Let us now turn to the study of the NNA predictions for the coefficients of the pertur- 
bative series for Rs in the minkowskian region. As in the case of Ds and D^-functions, 
this procedure gives the correct sign and order of magnitude at the 3- and 4-loop levels 
(see Table 4). Taking into account the large A^/-result for A5 (see Eq.(57)) we get the 
following NNA prediction for the 5-loop term in Rs 

- /'S (^ - ^C(3) - MC(5) - + 4C(3).= + i.') (107) 

which gives small and positive numbers 

fiNNA^^^ = 3) ~ 49 ; sT^{Nf = 4) ^ 39 ; s^^'^Nf = 5)^31 (108) 
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Following now the conservative pattern of fixing the uncertainty of the NNA approximation 
of dn-terms by an overall factor of 2, we present the NNA-inspired estimates of S4-term in 
the following form: S4 ~ 254^"^. Keeping this in mind, we arrive at the following numerical 
estimates of S4 for different numbers of Nf-. 

Si{Nf = 3) ~ 98 ; S4(A^/ = 4) ~ 78 ; s^{Nf = 5) ~ 62 (109) 

For Nf = 5 these estimates are quite in accord with the result of applying the [2/2] 
asymptotic Pade-approximation method, namely sf^^^{Nf = 5) ~ 67 p4|. However, for 
Nf = 3 the Fade estimate of Ref.[2^|, namely sf^^^{Nf = 3) ~ 251, is over 2.6 times 



larger than the related NNA inspired estimate of Eq. (|109|) . 

Note, that contrary to what was found in the application of the Fade-resummation 
method to the e^e~ annihilation i?-ratio [^, the variant of asymptotic Fade estimates 



used in Ref . p4| , which is performed in the minkowskian region directly, does not allow one 
to separate the effects of analytical continuation proportional to tt^. In the NNA approach 
the contributions of the vr^-effects leading in /3o are taken into account explicitly and can 
be resummed to all-orders of perturbation theory. We shall consider this technical problem 



in Sec. 3. 3. Another observation is that the estimates of Eq. ( |109| ) differ in both sign and 



order of magnitude from the ones obtained in Ref. p3|| using a variant of the effective- 



charges procedure [21, 63]. In Sec. 3. 6 we shall return to more detailed considerations 
of the problems related to the application of the effective charges method in the scalar 
channel. 

We conclude this subsection by demonstrating the behaviour of the MS- perturbative 
series for the physical quantity Rs of Eq.(71) in the cases of Nf = 3 and Nf = 5, which 
are relevant to the spectral function of the QCD sum rules and the hadronic decay width 
of the Higgs boson. Taking ag ~ 0.3 in the first case and ag ~ 0.114 in the latter one we 
have 

Nf = 3 : Rs ~ l + 5.667(^^) +31.864(^^)' + 89.156(^)' + 98(^^y (110) 
= 1 + 0.541 + 0.291 -1-0.078 + 0.008 

- 1 + 5.667(^^^ + 29.147(^^y + 41.758(^^^y + 62(^^y (111) 

= 1 + 0.206 + 0.0384 + 0.0021 + 0.00014 

One can see, that in both cases the perturbative series are rather well behaved and that 
the NNA-inspired estimates of 5-loop terms are over 10 times smaller than the 4-loop ones 
explicitly calculated in Ref.|I[] . In view of this we conclude that the 4-loop phenomenolog- 
ical studies, based on the 4-loop series of (110,111) are in rather good shape, and that the 
manifestation of the asymptotic growth of these perturbative series is postponed. Indeed, 
in accordance with the results of Tables 2,3 this feature can manifest itself starting from 
n = 7 loops. In its turn, this means that in the process of concrete phenomenological 
applications of the perturbative results for Rs in the energy region where ag < 0.3 one 
can restrict oneself to a consideration of the partial sums of the truncated perturbative 
series with n < 6 loops, estimating roughly the remaining perturbative uncertainty in the 
MS-scheme by the value of the smallest term taken into account. More detailed numerical 
studies are performed in Sec. 3.5. 
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3.2 Large Nj versus A^^ -theoretical motivation for NNA 



It is interesting to consider the motivation for the NNA approximation from the theoretical 
point of view. In the vector case the NNA term can be proved to have some very special 
properties by analyzing the operators that build the leading UV renormalon singularity 



65| . It is convenient to consider a 'planar approximation' where at each order in the 
Nf expansion only terms leading in Nc are retained. In this way one obtains an expansion 
in multinomials of Nf and Nc. So for the Adler-function coefficients in the vector case one 



can write, after extracting an overall factor of (3/4) [0, |6f 

dn = (t^N] + d^^~^^N]'^Nc + d^^-^^N]~^Nc^ + ... + S^'^NfNJ'-^ + rf^AT," , (112) 

so that the large- A^/ expansion runs from left-to-right, and the large- A^c from right-to- left. 
One can now formulate two versions of NNA. The standard one is derived by replacing 

Nf by (llA'c - 12/3o)/2 to arrive at 

dn = 4"^/3o" + 4""'^/3o"~'Ar, + 4"-2)/?o"-'A^c' + . . . + 4'^/5oiVc"-' + dfNJ^ . (113) 

One can, however , define a 'dual NNA' by replacing A'c by (12/?o + 2A'/)/ll, to obtain 
an expansion in (3q with different coefficients. The standard NNA is exact in the large- A^y 
limit, and the 'dual NNA' is exact in the large- Ac limit. Of course it is the standard NNA 
that is of practical use since we have all-orders large- A^/ results. If one extracts the NNA 
term d^^^ of the standard expansion and re-expands it one can obtain an expansion in 
Nf and A,, with coefficients d^^~^y By construction the leading- A^j term is reproduced, so 
d^n^ = d\^\ but the sub-leading terms will not be reproduced. Nonetheless by making use 
of the operator analysis of |Q one can show that [^] , 



4""'^^4""''[l + 0(l/n)], (114) 

so that for fixed-r and large orders of perturbation theory the re-expansion of the NNA 
term approximates the sub- leading in Nf terms with asymptotic accuracy 0(1/?t,). For 
the dual NNA term one can prove an exactly similar result, where sub- leading in A^ terms 
are reproduced to asymptotic accuracy 0(l/n) on re-expansion of the dual NNA term 
|65| . Such weak asymptotic results about the NNA terms will hold provided that in the 
large- A^/, and large- A^ limits the leading UV renormalon asymptotics is controlled by a 
single operator contribution. This is the case for the vector Adler function of Eq.(60). In 
planar approximation there are two relevant four-fermion operators, (9+ and (9_ of |]65| , 



but O- is scalar after Fierzing and decouples. These operators are defined as [^| 

where denotes the colour matrices. The remaining four-fermion operator (9+ gives 
the leading asymptotics in the large-A^c limit , and in the large-A^/ limit the operator 
corresponding to the single renormalon chains involved in NNA {Oi of [^) dominates the 
asymptotics IQ. This operator is defined as |^ 
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Let us see how these weak asymptotic results work by re-expanding the vector di,d2, 
expUcitly. We shall denote the dual NNA term by d^^"^*, the MS scheme with /i^ = is 
assumed, 



and 



d\ 

NNA 



NNA* 



-0.115iV/ + 0.655iVc 
-O.llSiV/ + 0.643iVc 
-0.119iV/ + 0.655iVc , 



(115) 



do = 



d 



NNA* 



0.086A^| - 1.40A^/A^e + 2.10A^c^ 
0.086A^| - 0.948A^/A^c + 2.61A^c 
0.069Ar| - {).7&2,NfNc + 2.10^^, 



(116) 



So we see that the weak asymptotic property holds, and the two versions of NNA give 
surprisingly good approximations for the sub-leading in Nf and Nc coefficients, given that 
this is only supposed to be an asymptotic result. 



Now let us repeat this analysis for the scalar £)5-function of Eq.(67). We find 

d2 



d. 



NNA* 
2 



= -l.dlNf + 17.19Nc 
= -l.dlNf + 10A9N^ 
= -3.13Nf + 17.19Nc, 



(117) 



and 



d 



ds = 0.93A^|-21.25A^/A^e + 72.08A^e^ 
0.93A^| - 10.22Ar/A^e + 28.10A^c^ 
2.38A^| - 26.21Ar^A^, + 72.08A^c^ 



NNA 
3 

NNA* 



'1181 



So we see that whilst the dual NNA term works reasonably well, the standard NNA in 
which we are interested yields the subleading in Nf coefficients with correct sign, but 
significantly reduced accuracy. For the coefficients of more physically-interesting quantity 
Rs of Eq.(71), which as we saw from the analysis of Sec. 2. 8, is perturbatively better 
behaved, we find 



S2 = 



^NNA 



^NNA* 
*2 



-1.36A^/ + 11.98A^c 
-1.36A^/ + 7.47A^c 
-2.18A^/ + 11.98A^c , 



(119) 



and 



S3 = 0.26A^| - 8.59A^/A^c + 18.24A^e 
0.26A^| - 2.85NfN^ + 7.S3N^^ 
OmN] - 6.63NfNc + ISMNc 



NNA 



NNA* 



(120) 
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We see that, again, the dual NNA term works quite well, but the standard version performs 
less satisfactorily. As in the vector case, this can be understood in terms of the operators 
involved 1^^. In the scalar case it is the four-fermion operator (9+ which previously 
dominated the large-iVc asymptotics which is vector after Fierzing and decouples. The 
remaining four-fermion operator (9_ will dominate the large- A''^ asymptotics, underwriting 
the success of the dual NNA. In the large- A^j limit, however, it turns out that (9_ and the 
one-chain operator in the scalar case are both involved in determining the asymptotics, 
and so standard NNA will not satisfy the weak asymptotic result that sub-leading in Nf 
terms are reproduced. Whilst the special property of NNA which holds for the vector 
case will not be true for the scalar, the numerical accuracy of the approximation is not 
so bad (see Table 4) despite the less satisfactory performance evident from the results of 
Eqs.(117)-(120). 



3.3 Analytic continuation of fractional powers of ag 

In this section we shall study the analytical continuation of the euclidean construct Ds{Q'^) 
introduced in (64). This is related to the quantity Rs by an analytical continuation to 
the minkowskian region. The continuation is essentially the same as that involved in the 
vector case for the analytical continuation of the Adler D-function to the QCD R-idXio. 
The effects of analytical continuation have been much studied H^-Q] in attempts to 



improve the convergence of the QCD perturbation series by resumming an infinite subset of 
analytical continuation terms at each order in perturbation theory. Such resummation may 
be accomplished conveniently by representing the continuation as a contour integral around 
a circle in the complex —Q^ plane (for one of the first discussion of this realization of the 



contour-improved technique see Ref.|^|). One can then perform the contour integration 
numerically , at some given order of perturbation theory. In the process one resums an 
infinite subset of potentially large analytical continuation terms involving powers of vr^, 
which arise in the running of the coupling around the circular contour. Such an expansion 
is termed "contour-improved" . For the case of a one-loop coupling an explicit closed-form 
result can be given for the contour integral. We would like to generalize these results to 
the present case where the mass anomalous dimension gives rise to fractional powers of 
as (recent analogous independent considerations were given in Ref . ||43|| ) . Using the NNA 
all-orders results for Ds and Rs we shall then perform various numerical studies on the 
performance of fixed-order perturbation theory, and its "contour-improved" version. 

The analytical continuation between D and R of (67) and (71) can be written in the 
form 

Rs{w') = ^ fje 3[m(eV)]' [l + Y.dn(^^^^^J^ , (121) 

involving a contour integral around a circle in the complex w"^ = —Q^ plane, as mentioned 
above. 

We can relate the running mass m{Q'^) to the renormalization scheme invariant mass 
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m as follows (see e.g. HI 



m{Q'^) = m exp 



MQ') -f^{x) ^ 70 ■ 



:i22] 



where the second term in the exponent is the commonly used normalization of the defini- 
tion of the invariant mass m. Since we shall be working within the NNA procedure, we 
shall set 7j = 0, (i > 0) and = 0(i > 0). In this approximation one has 

[m{Q^)f = m2(2/5o)'^»/^°(^^^y^°^* . (123) 

Inserting this expression for the running mass into (121) one arrives at 

n,i.,>.S.M-''-iiy[^ (124) 

The contour-improved expansion is obtained by performing the integration term-by-term. 
For a one-loop beta-function, appropriate for the NNA approximation, one has 

^s{e''^') = , , (125) 

[l + i/3oeas{w^)/-7T\ 

and so dn in the "contour-improved" NNA expansion will be multiplied by 



The function An{as) is given in closed form by 

An{as) = -T^ — (l + /^o«^)"^ (— ) sin [5„arctan (/^otts)] , (127) 



where (5„=1— n— 27o//3o. For 5„ — this reproduces the well-known factor (l/7r/3o)arctan(/5oQ;s 
obtained by resumming all the analytical continuation terms only involving f3o for the 
e'^e~ i?-ratio, while for n = the expansion of the r.h.s. of Eq.(127) to first order in 
as = 7r/(/3o/n(w^/A^)) coincides with the result previously obtained in [Q. Finally we can 
write down two expansions for R with NNA, 



or, alternatively, the "contour-improved" NNA expansion, 

270 //3o 



:i28i 



TT 



n>0 
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The A^'^^ for n > 1 are defined from (127) on setting 5„ = 1 — n. For n = one needs 
to be careful. The NNA terms are of the form /3o~^a*, with n = one needs to isolate the 
terms linear in 70 that arise on expanding (127) in powers of a^. One finds 

A^^^(a,) = 1 - ^ln(l + (3oW) - -^arctan(/5oa.) + • (130) 

Note that this term contains all contributions depending on the anomalous dimension 70. 
The remaining A^^^ for n > are precisely the same as the functions which arise in 
the case of the e~^e~ i?-ratio. The n = 1 case corresponds to 5„ — and so one has the 
well-known arctan alluded to earlier, 

A™(a,) = ^arctan(/5oas) • (131) 
TrPo 

For n > 1 the A^^^ are in fact simple rational functions of Os/tt. One has, for instance, 

aNNA/ \ {o^s/tt) 

A^'M = ((^) -^f^)1/(l + «)^ (132) 




3.4 Scheme dependence of NNA results 



Before proceeding to some numerical studies we need to confront one further important 
subtlety. As we have defined them the NNA expansions in (128) and (129) are scheme- 
dependent. Of course, we expect the partial sums to be scheme- dependent. The problem 
is that the all-orders sum for Rs will depend on our choice of renormalization scale. Since 
Rs is a physical quantity this is clearly undesirable. Following other similar analyses 
for the vector correlator [36, |6^, 69, |70[ we will use in the next section a Borel sum of 



the divergent series to define the all-orders sum, using regulation to cope with the IR 
renormalon contributions. The all-orders sum so defined can then be compared with 
fixed-order perturbation theory partial sums to obtain an estimate of the likely effect of 
uncalculated higher-order corrections. The problem is that the all-orders (Borel) sum of 
the series in (128) combined with the fractional power of (as/7r)^'^°''^° depends on the 
renormalization scale used for as- The difficulty is the fractional power of as involving 
1/13q. For illustrative purposes suppose that we used the so-called l^-scheme (see e.g. 



Ref . P9| ) corresponding to MS with jjr 



-5/3^2 



rather than /i 



w 



Writing and 



for the two scale choices, we have, assuming a one-loop beta-function. 



a 



MS 



270 //3o 



a 



y\ 27o//3o 



TT 



TT 



V 



1 + 

3 vr 



-270 //3o 



a 



v\ 270 /A) 

TT ) 



+ . . . 



(133) 
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Recalling that the NNA terms have the structure /9o~^«!, it is clear that only the terms 
linear in 70 in the expansion in the second line, appear in NNA. Thus, since not all terms 
involved in the change of scheme are retained, the resummed NNA expansions will be 
scheme-dependent. The resolution of this problem is to avoid powers of involving l//3o- 
This can be accomplished by identifying an effective charge R related to R by 



Rs = 37h^[2(3oR] 
R will have the perturbative expansion 



270 /A) 



R 



TT 



(134) 



(135) 



Since this only involves integer powers of all terms involved in a change of scheme at 
the one-loop level now contribute to the NNA result, and the resummed NNA expansion 
will be scheme- independent. 



Rs = 3m\2Po) 



270 //3o 



IT 



n>0 



TT 



270 //3o 



(136) 



is a reformulation of NNA for Rs in which the resummed series is scheme-independent. The 
unsatisfactory "scheme-dependent" version in (128) follows from it if only the terms linear 
in 7o are retained in expanding the series in sf^^^. Writing 3'^""'^^° = exp[(27o//3o)ln(5')] 
and expanding the exp to 0(70) one arrives at 



n>0 \ 



NNA ^s{w )\ 



TT 



1 J. 270, 



1 + E 



,NNA as{w )\ 



n>0 



TT 



(137) 



which relates the two versions of NNA for Rs- Using this result one can rewrite the 
reformulated expansion of (136) in terms of the sf^^, (128) is replaced by. 



Rs = 3m\2Po 



,27o//3o / asjw^y 



270 //9o 



TT 



exp 



E 

.n>0 



NNAf O^siw'^y 



IT 



:i38) 



We can immediately write a contour-improved version which replaces (129), 

270 //3o 

TT 



R, 



3m^(2/?o)^-/^°f^' 



exp 



aNNA/^ f„,,2\\ -r I sr^ ,NNA aNNA/^ /'„,,2\\ 



n>0 



(139) 



Using the Borel Sum to resum the series in the exponent in (138) and combining with the 
fractional power of as we will now obtain an all-orders result for Rs which is independent 
of renormalization scale, as required since Rs is a physical quantity. 
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3.5 Numerical studies on the convergence of the NNA results 



We shall now perform some numerical studies on the reformulated NNA expansions of 
(138) and (139). We shall consider the partial sums and Ry^ obtained by summing 
the series in the exponent in (138) up to and including the s^'^^ term, in the MS and 
V schemes, respectively. The prefactor 3m^(2/?o)^^'"'^° is set to unity. We shall also con- 
sider the analogous partial sums of the contour-improved expansion in (139), R^^^ and 
^(n)ci'^ In Table 5 we begin by displaying the s^^^ coefficients in the MS scheme, and the 
V"-scheme. We assume Nf = 5 active flavours of quarks. We shall use these two schemes 
to illustrate some of the scheme- dependence subtleties discussed in the last Section. In 
contrast to the "plateau of tranquility" evident in the limited growth of the MS scheme 
coefficients for n < 7 , which was alluded to in Sec. 2.9, we see that the corresponding 
V^-scheme coefficients grow much more rapidly. Alternating sign growth is evident even 
in low orders reflecting the asymptotic alternating factorial growth contributed by the 
leading UV renormalon singularity at 5 = — 1 contained in the Borel transform G^{6) 
in (44). In the MS scheme this behaviour is temporarily screened in low orders by the 
exp(55/3) factor in (53), which is absent in the V^-scheme. In Table 6 we show the par- 
tial sums for the choice of coupling = 0.114, (a^ = 0.12895) appropriate for Higgs 
width determination {Nf = 5 is assumed). As can be seen convergence is rapid for all 
the expansions. Of course, this is only temporary since the series is asymptotic, and for 
sufficiently large orders the alternating factorial growth of coefficients due to the leading 
ultra-violet renormalon will be evident. The \^-scheme leads to slightly faster convergence 
than MS. We show the s^'^^ coefficients in the MS and ^-schemes up to tt, = 10. Further 
note that for ra > 5 the partial sums in the two schemes are in complete agreement to the 
number of significant flgures quoted, emphasising the scheme- invar iance of the resummed 
expansions in (138) and (139). In Table 7 the partial sums for af^^ = 0.3, (a^ = 0.46736) 
appropriate for the strange quark mass determination are recorded {Nf = 3 is assumed). 
As one would anticipate the convergence is much less impressive. Further the analytic 
continuation terms are much more important at this larger value of the coupling, and the 
agreement of the results in the two schemes, and the apparent convergence is much more 
evident for the contour-improved expansion (139). This fact supports the application of 
the contour-improved NNLO expansions for the extraction of the s-quark mass value from 
the Cabibbo suppressed r-decay mode directly in the MS-scheme and within a re- 
alization of the effective charges approach To emphasise the scheme-dependence of 
the all-orders sum of the "conventional" NNA expansions in (128) and (129) we tabulate 
the corresponding partial sums in Table 8, for af ^ = 0.114. The partial sums in the two 
schemes are clearly converging towards two different results, 0.04134 in the MS scheme 
and 0.04039 in the l^-scheme. The difference ^0.001 is of order (as/vr)^, as one would 
anticipate from (133). 



We can compare the partial sums in Tables 6-8 with the all-orders Borel sum of the 
series based on the Borel transform in (53). IR renormalons require regulation since they 
contribute singularities on the positive real axis in the Borel plane. In common with 
similar numerical studies on the vector correlator p6|, |66|, |69|, [70[] we shall take a Cauchy 
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Principal Value (PV). In the V^-sclieme we then have 



E 

71>0 



d6 e-^"/(*"^) 



(140) 

where the UV and IR renormalon contributions G-{6) and G+{6) are given by (44) and 
(45) respectively. Writing the 'sin' as a sum of complex exponentials and using partial 
fractions the separate UV and IR renormalon contributions can be expressed |^6[ in terms 
of generalized exponential integral functions Ei{n, w), defined for Rew > by 



Ei(n, w) 



OO g 

dt 

1 



wt 



:ui) 



One finds 



- 

71 J 



A 



NNAf^y 







l + 4A™(ar) 



4 i-lY 

^^^[0+(l, n) + 0+(2, n) - 0_(1, n) - 0_(2, n)] 



+ 



3vr/?o 



7r/?o 



n>0 



f^-(l.l) 



where and (j)- are defined by [|6^ 



-(P,g) 



^0_(l,2) + ^0_(2,l)-^0_(2,2) 



(142) 



g-W(/3o"r)(_i)97r / g 



p-1 



(p-1)! 



TT 



P-1' 



(143) 



Including 1000 terms in the sum over 0+ and 0_ in (142) gives a result accurate to five 
significant figures. Exponentiating and evaluating R from (138) then yields the values 
indicated "P\/" in the last row of Tables 6 and 7. The all-orders Borel sum in (140) is 
scheme-dependent by itself, but on exponentiating it and combining with the factor of 
one obtains a scheme-independent result by construction. In the MS 
scheme the Borel transform in (140) has an extra factor e'^/^. At the smaller value of the 
coupling af^ = 0.114 in Table 6 one sees that for n > 4 the partial sums are all in good 
agreement with the exponentiated PV Borel result. The partial sums are in fact stable 
to four significant figures up to n^40 where violent oscillations due to the leading UV 
renormalon will be evident. The situation is somewhat less stable at the larger value of 
= 0.3 in Table 7, with oscillations due to the leading UV result clearly visible in the 
V^-scheme. For n > 4 the MS result is in reasonable agreement with the exponentiated 
PV Borel sum, oscillations become evident for n > 9, the ^-scheme results break down 
for n > 7. The contour-improved expansion results are stable but somewhat smaller 
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than those obtained with the standard perturbative expansion. On the evidence of these 
numerical comparisons one would anticipate that, even at the larger value of the coupling, 
fixed-order perturbation theory at NNLO (n = 3), the level of exact calculation at present, 
would give a reasonable approximation. Finally, in Table 8 we give in the last row the PV 
Borel sums in the MS and V schemes. We see that these values are in good agreement 
with the respective scheme- dependent "conventional NNA" partial sums for n > 4. 



3.6 Effective charges from the scalar correlator 

In Sec. 3.1 we discussed the problem of estimates of the higher-order QCD corrections to 
the spectral function of QCD sum rules in the scalar channel and to the decay width 
of a Standard Electroweak Model Higgs boson to quark-antiquark pairs within the NNA 
approach. In fact this problem was already analysed a few years ago [^| using a variant 



of the procedure developed in Refs.[^ which is based on application of the effective 
charges approach of Ref.|^. Usually this approach is applied to the renormalization 



scheme-dependent expansions of the quantities, which satisfy the renormalization group 
equations without anomalous dimension terms (see e.g. Eq.(66)). However, the quantities 
of Eqs.(67),(71),(91) we are interested in obey the renormalization group equations with 
anomalous mass dimension function (see e.g. Eq.(49)), which arises due to the factor 
of two powers of the running quark mass appearing in the Born approximation of their 
expansions. The appearance of the anomalous dimension function in the corresponding 
renormalization group equations reflects the scale-scheme dependence of the running quark 
mass and it generates the additional scheme-dependence of the perturbative series under 
investigation (in contrast to the familiar e~^e~ annihilation R-ratio, the scheme-dependence 
of Rs is starts from the as/vr-term). This additional scheme-dependence is analogous 
to the factorization scale-scheme dependence of the moments of deep-inelastic scattering 
structure functions (for some related discussions see Sec. 3. 4). In the process of the effective- 
charges motivated studies of Ref.|23| the careful treatment of this "factorization-like" scale- 



scheme ambiguity of the definition of running mass was overlooked. To overcome this 
shortcoming one should define the related effective charges either using the representation 



of Eq.(123) of Sec. 3. 3 p^, |63] or defining the logarithmic derivative of the quantities under 



consideration 



In this section we shall follow the latter prescription and consider the euclidean con- 
struct -R^i in (^) which is the unique combination of first and second derivatives of 
the high-energy scalar correlator that satisfies a renormalization-group equation of the 
form ( |66[ ) and gives Rd = tts/vr + 0(a;^). Thus it provides the simplest way of defining an 
effective coupling 

= Rd{Q ) = --—- — -—— = ^rn\ 144 

TT 2 dlogQ2 y ^ J 

in the scalar channel. This leads to an effective beta function 

~ rflog^B(Q^ /a^(g2)y+^ 
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whose coefficients /5„ are scheme-independent combinations of the coefficients (3n of the 
MS beta function and the coefficients r„ of the MS expansion of Rd- Clearly may differ 
from Pn only for n > 1. To simplify the presentation, we define Cn = Pn/Po- Then for 
n > 1 we have 

Pn Pn 



n 



/3o(r„ - 



(146) 



where Qn is determined by products of elements of {r^, Ck \ k < n} with weights summing 
to n. In particular 



^^2 
^4 



ri\C2- -CiTi 



2r{ + 3r2 

2 14 
ri I C3 - -Csri + -cira + — 
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-rir2 + 4r; 



+ 2 ( C2r2 - cirs + 5r2 



(147) 
(148) 



(149) 



to 5 loops (see |^T], ^ 



The estimation method based on this effective charge is f» VL^. It was motivated by 
the assumption that in electron-positron annihilation and deep-inelastic scattering sum 
rules such scheme-independent process-dependent effective P functions will behave in a 
way that is broadly similar to the MS P function, with /3„ having the same sign and 
magnitude as /3„ at 3 and 4 loops. The calculation of Ref. [1^ confirmed this assumption 
(for more detailed comparison of the behaviour of the effective charges /^-functions for 
these processes with the MS /^-function at the 4- loop level see Ref.p6|). 

Note however that the effective coupling R£,{Q'^) contains the spurious IR renormalon 
at 5 = 1, because it was constructed from the ffist and second derivatives of the correla- 
tor, with the ffist requiring UV subtraction for finite rr? jQ^ . The only way that massless 
perturbation theory can tell us that we did something that would be illegal in the mas- 
sive theory is to make the massless perturbation explode factorially, leaving an intrinsic 
ambiguity of order A^/Q^ at the point where the sign-constant asymptotic series becomes 
senseless. We emphasize that this A^/Q^ effect is profoundly perturbative and should not 
be confused with suggestions of A^/Q^ effects beyond those expected from the OPE ||72 . 
In our present case we know precisely why Ds{Q'^) is sick: the affiiction results from per- 
turbative malpractice, in failing to remove UV infinities at finite mass. If one follows the 

by taking the second derivative, then the spurious 5 



1 renormalon 



good practice of |142[ 
disappears. 

It is therefore instructive to compare the performance of the Vt estimator with its 
obvious alternative ^ Q4 based on the effective charge 



Rn{Q') = -- 



2 dlogQ^ 



n>0 



71+1 



TT 



(150) 



where r„ is obtained by replacing dn in (p6|)- (|100|) by dn from (|9TD . This effective charge 
is the unique choice formed from the second and third derivatives of the correlator and is 
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hence the simplest that is free from the spurious renormalon. To compute f2„ exactly one 
merely replaces r„ by in ( |147] )- ([149|) . 

In Tables 9 and 10 we compare the estimators Q and Q with exact results at 3 and 4 
loops. In Table 9 we also give the naive nonabelianizations of r„ and fin- Those for Tn 
and Q„ are in Table 10. Here 



,NNA 



^NNA 



-/?0 



(151) 
(152) 



are obtained by Nf —>■ Nf — 33/2 in (|101|) and its corresponding version for the effective 
charge RoiQ'^)- The NNA expressions for the estimators Q are 



or 



^NNA 



^NNA 



,NNA 



4r. 



NNA^NNA 
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n 



NNA 



r. 



NNA 



(153) 



3r2™r™-2(^rf^^y 



5 / nna\ ^ 



^A3A^ + 1a^ + ^A^ 



5 
3 



(154) 



while the ones for fl can be obtained by rf^^ to ff^^ in Eqs.(153) or A„ to A„ in 
Eqs.(154). For example, the analytical form of the NNA expressions for estimators and 
read 



^NNA 



■NNA 



4?[C(3)F.>.1^ 

/3o^(f[C(3)F + ^C(3) + ^ 



(155) 
(156) 



They result from neglect of Ck in (|149|) , since Ck = 0{Nj ^), at large Nj, while and 
are 0{Nf). 

In Table 9 we present results using the Q estimator. The comparison with exact 



results at 3 and 4 loops is not so good as in the QCD and QED studies of Refs. pT|, |6c 
correspondingly. At Nf = 5, for example, we have r2/fl2 ~ 0.702, while r^/Qs ^ 2.34, so 
the exact growth factor r^/r2 is 2.34/0.702 ^ 3.33 times that estimated by Q. Considering 
the NNA approximations for r„ and Qn we get the similar results. Thus effective-charge 
analysis, based on both and Vl^^^, give the results over factor 3 larger than the growth 



factor estimated by ( |151 
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At Nf = 5, this gives r2^^/r2 ~ 0.60 and r^^^/r^ ^ 0.59, with an actual growth factor 
r3/''^2 only 2% greater than that indicated at large-/3o- This is in accord with the global 
pattern of Table 4, which shows that NNA gives a reasonable account of growth from 3 to 
4 loops in all 3 quantities that were considered there. Now we turn attention to the final 3 
rows of Table 9. The values of ^4 in the final row are exact, since by definition they entail 
only input from lower orders of perturbation theory. The value of r4 is quite unknown: 
the Q method takes ^4 as its estimate. The smaller values of r^^^ are simply obtained by 
replacing Nf by Nj — 33/2 in the exact large- A^/ result for r^. The intermediate values of 
^24^^^ come from (|155|) . The A^/-independent ratio r^^^/il^^^ ~ 0.46 is a precise measure 



of the uncertainty of the fl estimator at large Nf. The A^/-dependent ratio Q^^^/Q4^ ?s 0.5 
is a precise measure of the uncertainty of NNA for Nf = 3, 4, 5. These two effects lead to 
a factor of 4 difference between the estimators r^^^ and fl^, as rival candidates for r4. 

Now we turn to Table 10, where the Q method fails to get the sign right for rs, at 4 
loops. Inspecting the final 3 rows, we see that the fl estimator is 10^ times its target at 
large Nf. On the other hand, the NNA estimate (152) is very successful. In view of the 
these failings of Q, we proceed only with the Q estimator. 

To convert a prediction of into one for the 5-loop term S4 in the physically relevant 
imaginary part, we use the known terms in the quartic 

r4 - 74 - 2/?oS4 = 37136.85285 - 78 1 0.21 6455 A^/ + 575.3282994A'J 

- 16.32062026A^J + 0.1444281007A^J (157) 

and parametrize errors in the Q estimator by Vn/fln = 1 + 5n- Then we obtain 

S4 = S4- 74/2/^0 (158) 

s^(Nf = 3) = 472 + 4574^4 (159) 

s^(Nf = 4) = 146 + 3529^4 (160) 

s^(Nf = 5) = -129 + 2615(54 (161) 

with central values that are small compared with a realistic estimate of the uncertainty, 
bearing in mind that 62 ~ —0.3 and ^3 ~ 1.3, at Nf = 5. 

As an indication of the problem at 3 and 4 loops, we give the 4-loop effective beta 
function in terms of dg = Rd at Nf = 3, 4, 5: 

f3{Nf = 3) = -2.2505, - 4.0005^ + 58.9205^ - 2148.5035^ (162) 
j3{Nf = 4:) = -2.0835, - 3.2085^ + 53.852a^- 1687.191fit (163) 
f3{Nf = 5) = -1.9175, -2.4175^ + 49.356a^- 1303.4905^ (164) 

with a sign change at 3 loops as a result of the appearance of a large and positive 3-loop 
coefficient. This pattern was already observed to occur in the effective beta-function for 
the minkowskian analog of Rd in Ref. 0, where it was considered as an indication of the 
existence of the spurious perturbative infrared fixed point. Indeed, this zero is compen- 
sated by the 4-loop terms, which remove the spurious fixed point (for the demonstration 
of the existence of a similar feature in the minkowskian region see Ref . ) . 
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They are however 0(50) times larger than those in the MS-scheme beta functions 





= 3) = 


-2.250a, 


- 4.000a^ - 


lo.oeoa^ 


- 47.228a^ 


(165) 




= 4) = 


-2.083a, 


- 3.208a^ - 


6.349a^ 


- 31.387a^ 


(166) 


f3{Nf 


= 5) = 


-1.917a, 


- 2.417a^ - 


2.827a^ 


- 18.852a^ 


(167) 



with a, = a,/7r. 

From this perspective, it is unsurprising that the VL estimator performs less reliably 



at 3 and 4 loops than in the cases considered in Refs.[^, |6J: at 3 loops it is bound to 
overestimate r2, because of the sign change of /?2; at 4 loops it is bound to underestimate 
r3, because [3^/ [3^^ ~ 50, whereas the procedure assumes that it is 0(1). 

We conclude that the application of Vt estimator in the scalar channel allows values 
of the 5-loop coefficient S4 between —10^ and +10^, if the accuracy is no better than 
^4 = ± 0.3, i.e. no better than at 3 loops. If it is no better than at 4 loops, the range 
widens further, by a factor of about 3. 

The appearance of large and negative estimates for S4 in Ref. ||23|, which contradict the 
results of application of the NNA procedure, described in detail in Sec. 3.1, is a reflection of 
the similar problem encountered in Ref. |2^ in a simplifled variant of the effective-charges 
approach in the same scalar channel. 



4 Conclusions 



In this paper we have extended the existing large- A^j analysis of the vector correlator [14] 
to the previously uninvestigated case of the scalar correlator. Because of the absence of 
the Ward identity Zi = Z2 present in the vector case, and the inevitable involvement of 
the quark mass anomalous dimension, the analysis amd combinatorics was considerably 
more complicated. An all-orders large- A^y result for the anomalous mass dimension 7^ was 
obtained in (39,40), and thanks to the remarkable identity (20) relating insertions into two- 
loop skeleton diagrams to the anomalous dimension, it was possible to obtain the all-orders 
large- A^/ result for the coefficient function of Rs in (51)-(53). As in the vector case the 
Borel transform of (44) and (45) contained UV and IR renormalons. A new feature of the 
scalar analysis was the presence of a leading IR renormalon at 5 = 1 in (45). This is not 
present in the physical quantity Rs thanks to the analytical continuation factor sin(7r5) / tt6 
in (53), and the leading singularity is the UV renormalon at 5 = — 1. It is, however, present 
in the quantity DsiQ"^) of (64), which would correspond to the obvious generalization of 
the vector Adler function , related by a singly subtracted dispersion relation to the scalar 
vacuum polarization 115. The presence of this leading IR renormalon is connected with 
the fact that the undetermined constant n5'(0) in 115 is infinite, a circumstance which did 
not occur in the vector case thanks to the Zi = Z2 Ward identity. A more satisfactory 
choice for the scalar Adler function is therefore the twice subtracted construct Ds defined 
in (63). An initial survey of the growth of the coefficients S'„ in (50) was given in Table 
2, where it is seen that in the MS scheme for n < 7 there is rather stable behaviour, 
with rapid growth corresponding to the leading UV renormalon evident for n > 7. The 
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perturbation series for Ds,Ds and its logarithmic derivative Rd are studied in Sec. 2.10. 
In Sec. 3 we moved to a study of so-called "naive non-abelianization" (NNA) , where Nf 
in the large- iV/ result is replaced by Nf — llNc/2 = — 6/3o, and the pieces of perturbative 
coefficients containing the leading power of /?o are resummed to all-orders. For the quanti- 
ties Ds,Ds and Rs the use of NNA to estimate the known three and four-loop coefficients 
was found to give the correct sign and order of magnitude. Moreover, the NNA-inspired 
estimates predict that at the five-loop level the numerical values of the corresponding co- 
efficients are positive and not very large. In the case of the imaginary part of the scalar 
correlator, related to the Higgs boson decay width to quark-antiquark pairs, and to the 
QCD Sum Rules spectral functions, we have sf^'^iNf = 3) ^ 49, s^^'^iNj = 4) ^ 39, 
s^'^^^Nf = 5) ~ 31. Using the conservative estimate of an overall factor of two for the 
uncertainty of the NNA estimates in the scalar channel (suggested by careful comparison 
of the three and four-loop NNA estimates to the results of explicit calculation, see Table 4) 
we conclude that the NNA-inspired estimate ~ 2^^^"^ gives us the following numbers 
s^{Nf = 3) ~ 98, S4{Nf = 4) ~ 78, s^{Nf = 5) ~ 62. For Nf = 5 our estimate is in good 
agreement with the result of previous studies, based on application of the [2/2] asymptotic 
Fade estimation technique [24] . However, for Nf = 3 the method of Ref . [^ agrees with 
our estimates only in sign and is over 2.6 times larger than the estimates proposed by 
us. In Sec. 3.2 we noted that for the vector correlator the NNA term does have some 
special properties, deriving from an analysis of the operators that build the leading UV 
renormalon singularity. On re-expansion it reproduces the sub-leading in Nf contributions 
to asymptotic accuracy 0(l/n) in n^^ order perturbation theory. A similar result holds 
for a "dual-NNA" , exact in the large- A^^^ limit. For the scalar case a corresponding weak 
asymptotic result was not expected to hold, but the "dual NNA" term should still provide 
a good approximation. In Sec. 3.3 we considered the analytical continuation from Ds 
to Rs ■ We recast the running mass m{Q'^) in terms of the RG-invariant mass m . The 
analytical continuation was similar to the much studied Euclidean to Minkowskian contin- 
uation for the e^e~ i?-ratio [29-40]. The presence of an anomalous dimension complicated 
the analysis slightly. We arrived at the "contour-improved" expansion for Rs in (129), in 
which a subset of analytical continuation terms involving vr^ are resummed to all-orders, 
at each order of the expansion. A further subtlety due to the presence of an as^'^°^^° term 
involving the anomalous dimension, was that the straightforward NNA expansions in (128) 
and (129) have all-orders sums which are RS-dependent. We showed how this could be 
remedied in Sec. 3.4 , and obtained the reformulated expansions in (138) and (139) whose 
all-orders sums were scheme-independent. In 3.5 we performed comparisons of fixed-order 
perturbation theory with the all-orders sum defined using a Cauchy principal value of the 
Borel sum, based on the Borel transform of (44), (45), (Tables 6-8). Two values of the 
coupling, = 0.114 and af^^ = 0.3, appropriate for the calculation of the Higgs decay 
width to a quark-antiquark pair, and for the strange quark mass determination from QCD 
Sum Rules, respectively, were considered. The MS scheme and l^-scheme were used to 
illustrate the scheme-dependence issue. Even at the larger value of the coupling it seemed 
that satisfactory accuracy was achieved at order n = 3, the highest order at which exact 
calculations are so far available. However, from the analysis of Ref. |]73| one can conclude 
that the uncertainties of the existing rris extractions from the QCD sum rules for the 
scalar correlator might be underestimated. In view of our analysis we conclude, that these 
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possible additional uncertainties are not coming from the uncalculated 5-loop perturbative 
QCD contributions, but are mainly related to the uncertainty of the experimental model 
for the spectral function of the QCD sum rules in the scalar channel. 

In overall conclusion we have provided a framework for extending the numerous pub- 
lished investigations of higher-order perturbative behaviour, analytical continuation, ap- 
proximate all-orders resummation, and estimates of higher-order corrections, for the vector 
correlator and its derived quantities, to the much less studied case of the scalar correlator. 
We hope that our results may be of use in further phenomenological investigations. 



5 Note Added in Proof 



After this work was submitted for publication we were informed that the 3-loop corrections 
to the correlator of scalar currents have been calculated up to (m^/(5^)^-terms [73|. The 
application of Pade resummation technique |^5|, |7^ allowed the authors of Ref.[^ to 
specify the exeat mass-dependence of the scalar-scalar correlator at 3-loops in a semi- 
analytical way. The analytical mass dependence of the 3-loop double-bubble contribution 
to the scalar correlator was studied in Ref . ||78[| . It should be stressed, however, that none 
of these calculations affect the results reported in this our work. 



Another, more closely related investigation, was performed in Ref. , where the effects 
of the UV-renormalons to the IR-safe Adler function of Eq.(63) of the scalar correlator 
were studied both in the 1-renormalon and 2-renormalon chain approximations. However, 
this analysis was based not on the exact calculations of the renormalon chain diagrams in 
the large A^j-limit, but on the latge iVc- analysis of the contributing 4-fermion operators. 



Moreover, in Ref. |^ the analysis of the IR renormalon contributions was not considered. 
In view of these differencies, it would be of interest to compare our results with those of 



Ref. 791 in more detail. 
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Table 1: Renormalons in (^) 



n 


jjMS 


n 


Hi-) 

n 


n 


i/(2) 

n 


H(+) 

n 


Hn 


2 


19.625 


-12.500 


2.594 


-48.000 


9.750 


1.043 


-27.488 


3 


11.093 


6.944 


5.870 


73.000 


-11.125 


-1.306 


84.476 


4 


-8.684 


-5.787 


12.276 


-173.778 


18.069 


2.289 


-155.615 


5 


-0.056 


5.787 


45.579 


556.056 


-35.785 


-4.764 


566.817 


6 


1.882 


-6.430 


207.742 


-2270.519 


83.144 


11.212 


-1972.969 


7 


-0.509 


7.655 


1177.991 


11416.893 


-223.934 


-29.266 


12348.830 


8 


-0.142 


-9.569 


7883.559 


-68593.216 


694.766 


83.893 


-59940.709 


9 


0.095 


12.404 


60911.356 


480286.469 


-2465.171 


-262.531 


538482.621 



Table 2: Renormalons in (pl|) 



n 


CMS 
'-'n 


^(0) 
n 


n 


5(1) 


^(2) 






2 


-1.667 


3.333 


1.097 


4.000 


-1.000 


-0.097 


5.667 


3 


0.972 


-0.512 


-0.576 


10.667 


-2.167 


-0.232 


8.152 


4 


0.834 


-7.880 


-0.999 


19.285 


-1.655 


-0.263 


9.323 


5 


0.055 


-13.817 


-2.495 


10.576 


9.338 


0.762 


4.419 


6 


-0.155 


9.160 


10.385 


-68.229 


46.309 


5.344 


2.813 


7 


-0.056 


105.021 


-17.734 


-249.544 


92.873 


15.175 


-54.266 


8 


0.006 


200.490 


183.066 


-237.016 


-36.445 


13.543 


123.645 


9 


0.008 


-255.017 


-1342.424 


1147.753 


-829.274 


-87.712 


-1366.667 


10 


0.001 


-2257.800 


9978.990 


5101.747 


-2297.027 


-516.006 


10009.905 


11 


0.000 


-4046.531 


-91303.959 


5275.323 


-176.520 


-1500.459 


-91752.148 



Table 3: Contributions to (||,|9|,|3|) 



n 




A„ 


A„ 


A„ 


2 


-1.66667 


7.33333 


7.33333 


5.33333 


3 


0.97222 


10.4696 


7.17968 


3.13622 


4 


0.83422 


32.6145 


8.48885 


11.6754 


5 


0.05531 


97.9534 


4.36402 


0.10978 


6 


-0.15502 


503.887 


2.96849 


112.074 


7 


-0.05583 


2194.28 


-54.2101 


-325.157 


8 


0.00581 


16465.8 


123.639 


3300.11 
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Table 4: Ratios of NNA estimates to exact results 





iV/ = 3 


Nf = A 


Nf = 5 




0.514 
0.507 
0.498 


0.496 
0.490 
0.484 


0.477 
0.472 

0.469 




0.354 
0.482 

0.764 


0.346 
0.565 

0.871 


0.339 
0.747 

1.081 



Table 5: s„ coefficients in the MS and ^-schemes 



n 


MS 


V 


1 


7.3333 


4 


2 


13.761 


-4.3408 


3 


31.184 


6.7614 


4 


30.727 


-58.738 


5 


40.061 


456.674 


6 


-1402.21 


-4161.77 


7 


6129.65 


47635.6 


8 


-129864.7 


-640093.5 


9 


1.8231 10^ 


9.8022 10^ 


10 


-3.2028 10^ 


-1.6900 10*^ 



Table 6: Partial sums of (137,138) for af^ = 0.114 



n 




•-MS 












1 





04099 





04210 





04049 





04153 


2 





04174 





04179 





04152 





04166 


3 





04180 





04181 





04180 





04174 


4 





04181 





04180 





04180 





04179 


5 





04181 





04181 





04180 





04181 


6 





04181 





04181 





04180 





04180 


7 





04181 





04181 





04181 





04181 


8 





04181 





04181 





04181 





04181 


9 





04181 





04181 





04181 





04181 


10 





04181 





04181 





04181 





04181 


PV 





04179 





04179 





04179 





04179 
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Table 7: Partial sums of (137,138) for af^ = 0.3 



n 




TZ7\ — 

-^MS 




Tzr\ — 

Ky 


-^MS 


Ky 


1 





24971 





33333 





21619 





25609 


2 





28934 





29778 





25056 





26237 


Q 
O 


n 


OUUO ( 


n 


OU ( uu 





26816 





27759 


4 





30162 





29310 





27509 





27676 


5 





30180 





31222 





27853 





27606 


6 





30109 





28234 





27906 





27800 


7 





30144 





34521 





27765 





26881 


8 





30061 





21538 





27433 





27246 


9 





30191 





76056 





26774 





27747 


10 





29934 





01704 





26192 





25527 


PV 





30138 





30138 





30138 





30138 



Table 8: Partial sums of the "conventional" NNA (128,129) for af^ = 0.114 



n 




•-MS 












1 





03977 





04159 





03939 





04110 


2 





04034 





04133 





04018 





04121 


3 





04039 





04135 





04034 





04133 


4 





04039 





04134 





04038 





04133 


5 





04039 





04134 





04039 





04134 


6 





04039 





04134 





04039 





04134 


7 





04039 





04134 





04039 





04134 


8 





04039 





04134 





04039 





04134 


9 





04039 





04134 





04039 





04134 


10 





04039 





04134 





04039 





04134 


PV 





04038 





04133 





04038 





04133 
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Table 9: Exact values of r„ and Qn, with NNA estimates of and fl. 





Nf = 3 


iV/ = 4 


Nf = 5 


n 

„NNA 
' 1 


10.17 

8.25 


9.56 

7.64 


8.94 
7.03 


r2 

„NNA 
'2 


90.78 
53.00 
121.44 
68.06 


77.13 
45.44 
106.02 
58.35 


64.06 
38.46 
91.28 
49.39 


^^3 

„NNA 
'3 

^^3 

r^NNA 


1091.55 
557.25 
620.75 
188.77 


822.19 
442.36 
424.80 
149.85 


585.53 
344.46 
250.41 
116.69 


^NNA 

j|nna 

^4 


5020.88 
11020.06 
20583.50 


3690.50 
8100.07 
14704.33 


2643.84 
5802.84 
10024.21 



Table 10: Exact values of r„ and fin, with NNA estimates of r„ and fl. 





AT; = 3 


Nf=4: 


Nf = 5 


ri 


7.92 


7.47 


7.03 


ttNNA 

'1 


6.00 


5.56 


5.11 


r2 


36.53 


29.94 


23.52 


T^NNA 
'2 


15.88 


13.61 


11.52 


^2 


76.75 


67.34 


58.25 


^NNA 
"2 


36.00 


30.86 


26.12 


rs 


233.13 


138.29 


55.79 


T7NNA 
^3 


199.49 


158.36 


123.31 




-145.01 


-183.52 


-219.13 


3 


-146.21 


-116.07 


-90.38 


tjtNNA 

' 1 


5.63 


4.14 


2.96 


T^NNA 


5921.10 


4352.18 


3117.87 


^4 


6619.34 


4648.70 


3131.62 
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